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Abstract. We study a class of quasi-cyclic LDPC codes. We provide precise con- 
ditions guaranteeing high girth in their Tanner graph. Experimentally, the codes we 
propose perform no worse than random LDPC codes with their same parameters, 
which is a significant achievement for algebraic codes. 
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1 Introduction 

The LDPC codes are codes that approach optimal decoding performances, with 
an acceptable decoding computational cost ((H |2l [3l). In this paper we present a 
class of quasi-cyclic LDPC codes and we show that we are able to guarantee some 
relevant properties of the codes. Experimentally, their decoding performance is 
comparable with the performance obtained by random LDPC codes. 

Traditionally, coding theory is divided into two main research areas: algebraic 
coding theory and probabilistic coding theory. The former (|4|) deals with codes 
endowed with a nice algebraic structure, which allows both to study their internal 
properties and to have efficient encoding-decoding techniques, the latter deals with 
convolutional codes, which are very difficult to study and randomly constructed, 
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but having superior decoding performance. However, the rediscovery of the LDPC 
codes by MacKay ([3]) triggered a radical change in coding theory: now we have 
linear block codes that may reach decoding performance close to the upper bound 
given by the Shannon limit (['51). Dozens of papers have appeared since MacKay's 
paper, some of them trying to endow some structure (either algebraic or geomet- 
rical) on LDPC codes. However this has seldom been successful, because the 
structure brings a regularity in the parity-check matrix of the code, which natu- 
rally pushes towards the creation of many dangerous small cycles in their Tanner 
graph. It the object of this paper to propose a family of LDPC codes, possessing an 
algebraic structure but not suffering from the performance limitations common to 
other similar families. The family of quasi-cyclic LDPC codes are of great interest 
for the possibility of exploiting the structure of the parity check matrix to achieve 
very fast and efficient encoding and decoding. Unfortunately the BER/SNR per- 
formance of this class of codes is known to be worst than random generated LDPC 
codes in particular for medium/long length. 

The object of this contribution is to study the family of quasi-cyclic LDPC 
codes, and provide precise conditions guaranteeing high girth in the Tanner graph. 
Although several families of quasi-cyclic LDPC codes have been proposed, no 
general study on their girth properties have been published. Previous researches 
presented in the literature focuses on studies of the properties of particular classes 
of quasi-cyclic LDPC codes constructed from circulant matrices obtained from a 
monomial. It is our purpose to fill this gap, by formally classify all cases when 
cycles of length less than 10 may arise in the general case. Therefore, in this 
contribution matrices obtained from polynomial are also considered. The study is 
restricted to polynomials composed of two or less monomials, the reason for such 
limitation is the fact that circulant obtained from polynomial composed by three or 
more monomials internal cycles of length h 6 always exist. Hence such polynomial 
are of no interest if codes with higher girth is wanted. 

From the classification obtained, it is obvious how to identify necessary and 
sufficient conditions for any quasi-cyclic LDPC code to have girth at least 10. Var- 
ious constructions are presented that perform no worse than random LDPC codes 
with the same parameters, which is a significant achievement for algebraic codes. 

The remainder of this paper is structured as follows: 

• Section [2] provides notations, recall some relevant well-known facts and 
prove some simple preliminary statements; 

• Section [3] deals with theclassification of the cases when cycles up to length 
8 may arise for a rather general code family; 

• Section [4] improves results from Section[3]for the generic quasi-cyclic case; 

• In section \5\ the existence of short cycles is linked with condition on the 
polynopmial representation of the circulant matrices; 
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• In sections |7J [8] and [9] varius subclass of quasi-cyclic codes are studied in 
details and detailed conditions to avoid short cycles are given. The perfor- 
mances of such codes is compared with other codes, 

• Finnaly in section [10] comments, conclusions, and outline some further re- 
search are presented. 



2 Preliminaries and notation 

In this section some known facts are recalled, some notation are given and some 
simple statements that will be useful later on are proved. 



2.1 LDPC codes and Tanner graphs 

The parity-check matrix H = (hij) of any binary [n, k, d] linear code C may be 
represented by a graph T, known as the Tanner graph (1013). The Tanner graph is 
formed by two types of nodes: the "bit nodes" and the "check nodes". Bit nodes 
correspond to matrix columns and check nodes correspond to matrix rows, so that 
there are r = n — k check nodes and n bit nodes. We connect the check node i 
to the bit node j if and only if the entry h^j = 1. There is no edges connecting 
two check nodes or two bit nodes (this kind of graph is called a bipartite graph). 
In other words, H is the adjacency matrix of F. 

Example 2.1. An example of a binary LDPC code and relative Tanner graph can 
be seen in Figure\l\ 
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Variable Nodes 



101 100000 
101010000 
000000111 
0010100 01 








Figurd P. Parity dheek mafifixil and the associated Tanner Graph, the presence of 
a cycle is highlighted 



Now we introduce LDPC codes - Low-Density Parity-Check codes - a class of 
linear error correcting codes. Historically, these codes were discovered by Gallager 
in 1963 in his PhD thesis HJ. These codes were largely ignored, because of some 
implementation issues. In the 1990's they were rediscovered by MacKay Q and 
now the research continues vigorously, with dozens of papers published every year 
(HQ). 

Definition 2.2. An LDPC code is a linear block code for which the parity check 
matrix has a low density of non-zero entries. 
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A (c, s)-regular LDPC code is a linear code whose parity check matrix H 
contains exactly c ones per column and s ones per row. 

We do not specify what we mean by low density because it depends on the 
context. For example, for a typical (3, 6)-regular binary code (rate 1/2) of block 
length n, there are only three ones in each column of H and so the fraction of ones 
in this matrix is 6/n. 

The decoding algorithm for these codes is usually called the "sum-product al- 
gorithm". We summarize some properties of these codes: 

• the LDPC codes have excellent decoding performance, near to the channel 
capacity (EDM), 

• the sum-product algorithm is based on the probabilities received from the 
channel and it may be idealized as a belief -propagation algorithm (see [ 10]), 
with information being passed and updated by the bit nodes to the check 
nodes and vice- versa, in a loop; 

• the information propagation is heavily hindered by the presence of small 
cycles in the Tanner graph ( fTTTl ): 

• the best LDPC codes are irregular and are created by some random-walk 
optimization algorithm (Bll9l fT2l[T3l[T4l ): 

• there is no known algebraic class of LDPC codes that has performance com- 
parable to the best known LDPC codes. 

There are two serious issues for a code not possessing an algebraic structure: 

• the encoding process is computationally expensive, 

• it is very difficult to study its properties. 

There are many family of LDPC codes that have been proposed endowed with 
an algebraic (or geometrical) structure, but none of them has clearly shown, at 
present, a decoding performance comparable with the random LDPC codes (see 

US M M M M M ETJ El ) . 
2.2 Girth 

Definition 2.3. In a graph, a cycle is a path that starts from a vertex v and ends in 
v. The girth of a graph is the smallest of its cycles. 

Obviously the girth of a bipartite graph is always even. The girth is considered 
one of the important parameters of a LDPC code, it is commonly accepted that the 
presence of short cycles in the graph is one of the main parameters affecting the 
coding gain achievable by the LDPC code l23l . The dependency of the perfor- 
mances of a LDPC code on its girth distribution, in particular when small cycles 
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have been avoided, is still under debate since mathematical proof has not yet been 
obtained. In contrast with the deteriorating effect of cycles on the performance 
of the LDPC codes, Etzion et al. ll24l l25l proved how cycle-free graphs cannot 
support good codes. Still, simulations and applications have shown that the belief 
propagation algorithm is generally very effective, even in the presence of cycles in 
the graph j26l l5l. Nevertheless it is commonly accepted that the presence of short 
cycles in the graph is one of the main causes of reducing the coding gain achieved 
by the LDPC code [23 ], and so the girth is considered one of the significant param- 
eters of a code. 

For Tanner graphs of (c, s)-regular LDPC codes, it is possible to give upper 
bounds on the girth, see ll22l . 

Theorem 2.4. Consider a (c, s)-regular [N, K, d] LDPC code. Let R = N — K 
and a = (c — l)(s — 1). If the girth g = 2 mod (4), then g < 41og a r + 2, 
otherwise g < 4 log^ r + 4. 

When r = 404 and a = (3 — 1)(6 — 1) = 10 (as for our simulation), we get in 
the worst case g < 41og 10 (404) + 4 = 14, but in practice it is very difficult to 
find such codes with g > 10, so that ensuring g > 8 is already interesting. 

One of the most promising families of LDPC codes with a nice structure was 
proposed by Rosenthal and Vontobel ( |[22l ). These are based on Tanner graphs 
built starting from Ramanujan graphs and hence are guaranteed to have a very high 
girth. Unfortunately, their decoding performance have been questioned dl27l0 and 
it is not evident how an efficient encoding could be implemented. 

There is a family of LDPC codes, which has been proposed by Fossorier ([ 28 ]), 
which is particularly interesting for us, because they are quasi-cyclic and so their 
structure is quite similar to ours. With Fossorier's codes, it is easy to get a girth 
as high as 8 or 10. However, the construction by Fossorier does not provide codes 
whose Tanner graph has a girth higher than 12, as shown by Fossorier himself in 
the same paper. 

Another interesting construction has been provided in lf2Tl[29l . They can get 
very high girth, but there are some open problems, in particular on how to get an 
efficient encoding. 

To simplify the search for cycle in a Tanner graph representing a H parity check 
matrix, a novel and convenient definition of cycles of length l c in an arbitrary binary 
matrix is given here. 

Definition 2.5. Let s,N,M be natural numbers with s > 2 and N, M > 3, and 

define l c = 2s. Let B be any N x M matrix over Z 2- A sub-set V ofl c entries of 
B is called linked if the entries lay in s columns and s rows. 

Definition 2.6. Let B be any N x M matrix over Z 2- A sub-set V of l c entries 
(l c = 2s) of B is called a 2s-cycle if: 

• V is linked, and 
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for any r such that 2 < r < s there is no linked sub-set W C V of2r entries. 



In order to clarify the meaning of the definitions two examples are shown in 
Figure l2~2l 



8 entries linked set that 
does not form a 8-cycle 



i entries linked set that 
does form a 8-cycle 



Figure 2: Example of linked sets with s = 4, t = 8 

Matrix (a) contains 8 entries, they form a linked set since they lie in 4 rows 
and 4 columns but they do not form a 8-cycle since they can be grouped in two 
smaller linked sets of 4 elements each. Matrix (b) represents a linked set of entries 
that form a 8-cycle, since the 8 entries lie in 4 rows and 4 columns but there are no 
smaller linked sets. 

Note that if a matrix column contains a point of a 2s-cycle V then it contains 
exactly two points of V. The same is true for the rows. Moreover a linked sub-set 
V of 2s entries either is a 2s-cycle or it contains at least a 2r-cycle with 2 < r < s. 



2.3 Circulant matrices 

Binary circulant matrices are important as they form the "bricks" with which parity- 
check matrices for quasi-cyclic LDPC codes are "built". 

Definition 2.7. Let m > 6. Let C be anmxm matrix over F2. C is circulant if its 
rows are obtained by successive shifts (to the right). The matrix C is weight-l if the 
weight of any row is I. In case of circulant ofweight-2, that are used extensively in 
this work, the polynomial representation of the first row, p(x) G F2[x], is called the 
polynomial of C. Let p(x) = x a + x b , with a < b. s(p) = min(6 — a, a + m — b) 
is called the separation of p (or ofC). 

Consider, for example, the following weight-2 circulant matrix. 



/ 1 


1 








\ 





1 


1 














1 


1 














1 


1 


V 1 











1 / 
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Its polynomial is p = x + 1, with parameters m = 5, b = 1, a = 0, = 1. 

Note, similar matrices can be described as an superimposition of two permu- 
tation matrices used by other authors (e.g. QUI ), in such case the exponents are 
related to the power of the single permutation matrices, and the analysis presented 
in this chapter can be rewritten with such notation. 

When the natural number m is greater or equal to 3 the equations: 

(2) a = b mod (m), p(x) = q(x) mod (x m + 1) , 
will be abbreviated with, respectively, : 

(3) a = b, p(x) = q(x) , 

where the polynomial congruence is in Z 2 [x] . 

Let p = x a + x b , sometimes some statements where the role of a and b may 
be exchanged are needed. To provide a concise formulation in these cases, the 
notation e(p) in congruences modulo m is introduced. Let / be any function / : 
N^Z, then 

• /(e(p)) = I, means "/(a) = I or f(b) = I", 

• /(e(p)) # I, means "/(a) ^ I and f(b) # I". 

In the case of a circulant matrix of weight-1, e(p) refers to the exponent of the 
monomial. 

This notation is extended to the case when more polynomials, p\ . . .p s , are 
involved in a function, as follow. Let p\ = x ai + x bl , . . . , p s = x as + x bs . Let / 
be any function / : N s 1— > Z. Then 

• f(e(pi), . . . , e(p s )) = U means that there is a combination (z\, . . . , z s ) 
where %i G {oj, b{\ such that f(z\, ■ ■ ■ ,z s ) = I and 

• f(e(pi), . . . , e(p s )) ^ ^> means that for every possible combinations (z%, . . . , z s ) 
where z { G {a i; bi} then /(zi, . . . , z s ) ^ Z. 

Remark 2.8. To avoid ambiguity, statements of kind: f(e(pi), ■ ■ ■ e(p r )) = /(e(pr+l)> • • • > e (Ps)) or 
/(e(pi), . . . e(p r )) # /(e(p r+ i), . . . , e(p a )) . 
will never be used. 

Some simple facts on weight-2 circulant matrices are collected here. These 
follow directly from the circularity of the matrix. 

Proposition 2.9. Let C = {cij} be an (m x m) weight-2 circulant matrix over Z 
2[x] with polynomial p(x). Then: 

1. Cij = 1 if and only if j =i + e(p), 

2- c X)V = Ct, y = 1 (with x 7^ t) if and only ift — x = ±s(p), 
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3- c XtV = c X;Z = 1 (with y z) if and only if z — y = ±s(p), 

4- c XtV = ct <y = 1 and c X;Z = c WtZ = 1 (with y ^ z, x ^ t, x ^ w) if and only 
if x — t = ±s(p) and w — x = x — t = ±s(p). 

Lemma 2.10. Let p and q be two polynomials in Z 2 [x] with degree at most m—1. 
Then 

s(p) = s(q) ^ s(p) = ±s(q) 

Proof. Note that s(p) and s(q) wee not greater than m/2 and positive. □ 

Definition 2.11. Let m > 3. Let p and q be two polynomials in Z 2 [x] with degree 
at most m — 1. We say that p is a shift of q if there is < fx < m — 1 s.t. 

p = x^q . 

In this case we write p ~ q. 

Observe that ~ is an equivalence relation in the set formed by all polynomials 
over Z 2 with degree less than m. 

Lemma 2.12. Let p and q be two polynomials in Z 2 [x] with degree at most m — 1. 
Then 

s(p) = s(q) 4^ s(p) = ±s(q) 

Proof. Note that s(p) and s(q) wee not greater than m/2 and positive. □ 

There is a link between the separation of a polynomial and its roots. 

Fact 2.13. Let m > 3. Let p and q be two weight-2 polynomials in Z 2[x] with 
degree at most m—1. Then 

p ~ q <^ s(p) = s(q) . 

Moreover, if p and q are both maximal or if they are both minimal, then the 
non-zero roots of p and q are the same (and with the same multiplicity) if and only 
ifp ~ q. 

Proof. Let p = x a p + x bp and q = x aq + x b i, with a p < b v and a q < b q . 

p ~ q s (p) = s(q). 
If p ~ q, thenp = x^q, for some < /j, < m — 1. That is, p+x^q = X(x m + 1), for 
some A <G !*2[x]. But deg(p + x^q) < 2m — 2, since deg(p) <m — l,dq<m — l 
and /j, < m — 1. Also, either d (X(x m + 1)) > m or A = 0. 
There are three cases: 
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1. deg(x M g) < m — 1. 

Then A = and p = x^q, which means 

b p - a p = ±s(p), s(p) = s{x^q) = s{x aq+ ^ + x bq+ ^) , 

s{x aq+ » + x bq+ ») = ±{(b q + /J,)- (a q + M )) = ±(b q - a q ) = ±s(q) , 
so that s(p) = ±s(g) and hence s(p) = s(q) (Lemma [2.12l ). 

2. deg(x ti q) > m and we have x^q = x aq+fl + x bq+fJ- , with a q + fi < m — 1 
and b q + \i > m. 

Then x b "+^ = x b i+^- m ( x m + 1) + x b i+v- m , so that A = x b " + ^- m and 
either a p = b q + fi — m, b p = a q + yu, or b p = b q + fj, — m, a p = a 9 + /U, 
which implies s(p) = ± ((a g + /i) — (b q + fi — m)) = ±(a g — b q + m) = 
= ±(6 g — a g ) = ±s(g) and hence s(p) = (Lemma [2.12l ). 

3. deg(x M g) > m and we have x^q = x aq+fl + x bq+fJ- , with a q + > m and 
^ + ^ > 

This case is the same as case 1), with the role of p and q exchanged. Since 
~ is an equivalence relation, we do not have to deal with it. 

s(p) = s(q) =4> p ~ q. 
If s(p) = s(q), there are four cases: 

• s (p) = bp — a p and s(q) = b q — a q . Then b p — a p = b q — a q . We may assume 
b q > bp, so that b q — b p = a q — a p , i. e. q = x bq ~ bp p. 

• s (p) = m — a p + b p and s(q) = m — a q + b q . Again b p — a p = b q — a q and 
so we may argument as before. 

• s (p) = bp — a p and s(q) = m — b q + a q . It is enough to take fi = b p — a q = 

a p + m — b q . x^{x aq + x bq ) = x bp + x a p +m = x b p + x ap . 

• s(p) = m — bp + a p and s(q) = b q — a q . Same argument. 

We now suppose both p and q minimal and we want to show that p ~ q if and 
only if they have the same non-zero roots (with the same multiplicity). The case 
when they are both maximal is analogous and will not be shown. 
It is obvious that two polynomials p and q have the same non-zero roots with the 
same multiplicity if and only if p = x l q, for some i. Assuming both p and q 
minimal, we have p ~ q 43> p = x^q, and so our desired logical equivalence 
follows. □ 

3 Cycle configurations for generic matrices 

In this section some notations, facts and lemmas useful to identify which cycles can 
exist in a given matrix are presented. A rather general class of matrices is studied. 
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The general results obtained here will be specialized to the quasi-cyclic case in 
following subsections. 

For the remainder of this chapter, if not differently specified, m, a, are nat- 
ural numbers with m > 3, a, j3, 7 > 1, 

Definition 3.1. Given a matrix B over Z 2, B is said to be in M. m , a ,p,i if it may 
be written as 



B 



A 



1,1 



A 



l,a/3 



A 



07,1 



A 



ccy,ce(3 



where the Aij's are binary square matrices of dimension m. This decomposition 
is referred to as the standard decomposition of B in At m a / 9 -y. Any matrix Aij 
is called a decomposition sub-matrixes.). For any i in {1, ... , 07}, that the set 
{Aij I 1 <j < a/3} is called a decomposition row(d.r) of B. Similarly, for 
any j in {1, . . . , a(3}, the set {Aij | 1 < i < 07} is called a decomposition 
columnfj.c.j of B. 

The d.s.'s are defined unambiguously and the uniqueness of the standard 

decomposition in M. m ,a,p,~f is obvious. The term "standard decomposition" will 
be used rather than "standard decomposition in M. m: a,p,y'- If 5 G M- m ,a,p,^ B 
can also be viewed as: 
(4) 



B 



L 



1,1 



L 



7,1 



L 



J 4(r-l)a+l,(s-l)a+l 



(r— l)a+l,sa 



where L r s is a square matrix of sub-matrices, with dimension am, (1 < r < 7, 

1 < s < /?). 

Remark 3.2. T/ 1 // G .M,™ Q 3 7 has full rank and (3 > 7, ;Y represent a binary 
linear code with dimension (/3 — 7)7710 arac? length f3ma. The information rate is 



K 

N 



((3 — j)ma (3 — 7 



(3mct 



(3 



If H G A4 m Q /3 7 w not full rank the rate of the code is lower and the value pre- 
sented above can be considered as designed rate. Note that in some cases adding 
redundant rows, hence not having full, can be used as a method to improve perfor- 
mances tanks to the extra checks that a codeword has to satisfy. 

It is clear that any rate can be achieved with this code construction by choosing (3 
and 7 appropriately. 

Definition 3.3. Let B € M m ,a,p,~f an d let {A,j}i<i<a7,i<j<a/?/ orm its standard 
decomposition. Define the sets of indexes Z = {ii, ■ ■ ■ ,i r } C {1, . . . , 07} and 
J = {ji, . . . ,j s } C {!,..., a[3}. With the notation Bj^j the sub-matrix of B 
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defined by the d.r.'s in 2 and the d.c.'s in J is denoted. The sub-matrix Bj.j is said 
to be of type (r, s) and Bj t j is called a decomposition minorfii.m.j. Given two 
d.m.'s Bj,j and Cx',j', they are considered equivalent if it is possible to obtain 
one from the other by d.r. permutations or by d.c. permutations or by both. An 
equivalence class is called a configuration of type (r, s). 

Foe example the following two matrices are decomposition minors of matrix B 
in @] and are equivalent configurations since D' can be obtained from D by switch- 
ing the first and second rows and then the first and second columns. 





Lip 


£1,4 


£1,6 




£2,4 


£2,2 


£2,6 


D = 


£2,2 


£2,4 


-^2,6 


, D' = 


£1,4 


£1,2 






£4,2 


£4,4 


£4,6 




£4,4 


£4,2 


£4,6 



Note that the relation defined on d.m.'s is actually an equivalence relation, so that 
"an equivalence class" makes sense. 

The following lemma will be useful later on. 

Lemma 3.4. Let I be the m x m identity matrix over Z 2- Let B € «M m a p 7 
and let {Aij } form its standard decomposition. Suppose that there are i and j s. t. 
Aij = I. Let B X)V be an entry of B included in A{j. Then 



B 



x,y 



1 



x = y 



Proof. It is know that x = x' + (i — l)m and y = y' + (j — l)m, with 1 < x' < m 
and 1 < y' < m. The pair (x', y') represents the components inside Aij. But A^j 
is the identity, so that 



B 



1 44> x = y' 44> x = y 



□ 



Using this notation and the definitions of cycles on a matrix given previously 12. 6[ 
the following lemmas are obvious. 

Lemma 3.5. Let Bjj and Cp t ji be equivalent d.m.'s of a matrix B £ A^ mjCfij a i7 . 
Then Bjj (strictly) contains a 2s-cycle if and only ifCj t j does. 

Lemma 13.51 allows us to talk about "configuration of cycles"(see Def. I3.3I ). 
meaning equivalent decomposition minors that contain cycles of the same type. 
All the possible (d.m.) configurations with 2s-cycles are classified next. 

Lemma 3.6. With the notation introduced above, let B € M.m,a,p,T Then the 
only configurations of B that may (strictly) contain a 2s-cycle are oftyp^ 



\A 



'For brevity any configuration that is the transpose of another is omitted. 
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I A ■ ■ A ■ • I 

\-"l,3l"-l,32\ , 



(1,*) 

(2,2) 
(2,s) 



A • • ■ 


4- ■ 1 


4 ■ • 

"■11,31 


A ■ ■ 

■"■11,32 




A ■ ■ 

■"12,32 




... a 


^2 Jl 


... a 



3s 
12, 3 s 



(S- 1,8-1) 



1, S ) 



A 

A 
A 



-1,31 



11,31 



A 

A 
A 



n,3s-i 



-1,3s 



n,3s 



A 



A 



is-ljs 



^i,3i 



A 



H,3s 



A 



is, 31 



A 



Is ,3s 



Proof. Since a 2s-cycle V needs s matrix rows, it needs at most s d.r's, and anal- 
ogously for the d.c.'s. So configuration (s, s) is the largest that can occur. Since 
matrix rows (and matrix columns) can be grouped in d.r.'s (d.c.'s), any d.m. of 
(s, s) is possible. This proves the claim. □ 

To any d.m. configuration, one or more cycle configurations may be associated. 
To proceed it is necessary to characterize the cycle configurations arising from the 
previous lemma. In order to do so, a convenient notation for a cycle configuration 
is presented. Any d.m. configuration in the statement of Lemma 13.61 is a sub- 
configuration of an (s, s) configuration, which can generically be represented as 

T\ : i ■ ■ ■ Ti >s 
T St \ ■ ■ ■ T s ^ s 

where any Tjj is an A^ k, for some h and k. A numerical representation for cycles 
is adopted as follows. Let V be any 2s-cycle contained in the (s, s) configuration 
as above. For any matrix Tjj With t^j > the number of points of V contained in 
Tjj is denoted. With d.r. and d.c. permutations, it can be supposes that 

(5) ti,i > Uj, 1 < i,j < s, t 1)2 > £1,3 > > h, s , t 2) i > h t i > > t S; i . 
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A cycle presentation fulfilling conditions ([5]> will be called a (l)-presentation. 

It is obvious that for a 2s-cycle it must be: 



(6) 



l<i ,j<s 



2s. 



For example, a 4-cycle V, a 6-cycles F, and a 8-cycle V may be represented 



as 



2 2 




V 



2 

1 1 

1 1 





and all these representations are (l)-presentations. To clarify the implications of 
the representation consider the following cycle configurations: 



W 



w 



1 
1 




they do not form valid ( ^-representations since in W ti t \ is not the maximum of 
the elements and in W the elements of the first row are not ordered. 

Note that cycle V allows another (l)-presentation that is a column permutation 
of it: 

2 2 
V= 2 


For convenience, in a cycle presentation d.r.'s and d.c.'s containing only zeros 
will be dropped. The previous presentations may be written as follows 



V 



2 2 



V 



2 2 
2 



V 



2 2 
1 1 
1 1 



The notation presented above leaves some ambiguity. In fact, for example, 
ti i = 2 does not specify whether the two points in Ti j lie in the same cycle 
column or cycle row or in neither. 

Remark 3.7 (Transpose). Let D be a d.m. of a matrix B € M. mi a,p,j and D T 
its transpose. It is clear that D contains a 2s-cycle if and only if D T contains a 
Is-cycle. In the general case, it is not possible to obtain D from D T by d.c. or d.r. 
operations, so D and D T are not necessarily equivalent according to the definition 
given. However, if all cycle configurations associated to a given configuration 
are classified, then, automatically, all cycle configurations for its transpose are 
obtained (by transposing all of its cycle configurations). 

Two lemmas are provided next, these are applied in the analysis of the 2s-cycle 
configurations arising in Lemma [331 
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Definition 3.8. The row weight of a d.r. in a cycle configuration is defined as the 
sum of the tij 's in the d.r, and similarly for the d.c.'s. 

Lemma 3.9. In any cycle configuration, column weights and row weights are even. 

Proof. By transposing, the statement for d.c.'s is true if and only if it is true for 
d.r.'s. It is here shown for d.r.'s. Given a d.r., if the relevant tjj's sum to an odd 
number, then one point in not in a cycle row with another. But, by definition, any 
point in a 2s-cycle shares one cycle row with one and only one other point in the 
cycle, hence there cannot be d.r. with odd weight. □ 



Lemma [3791 will be applied many times. To simplify the relevant notation, the 
phrase "Lemma l3.9l r-1" will be used meaning "Lemma l3.9l applied to the first d.r.". 
Similarly with notations like "r-2", "r-3". The same notation is used the columns 
(e.g. "c-2"). 

Lemma 3.10 (Isolation). In any d.m. configuration D containing a 2s-cycle, the 
following situation cannot occur with 2 < r < s: 



D 



Di 
D 2 



where: 



Ti,i ••• Ti, 

1/ 'y ^ * * * 1/ 'y . 



D? 



T, 



r+l.r+l ' ' 



■,r+l 



1,8 



Proof. It is known, by definition 12.61 how in any row (and column) of a cycle V 
there must be exactly two cycle points. If one cycle point lies in a d.r. of the D\ 
then also the second cycle point of that d.r. must lie in D\. Hence, for every d.r. in 
D\ there are two cycle points, the same is true for the d.c.'s. So in D\ there are 2r 
points that lie in r d.r.'s and r d.c.'s , hence they are linked, but r < s and so this 
contradicts the definition of cycle given (Definitions 12.61) . □ 



To ease the reading in situations where lemma l3~T0l is not satisfy for a certain 
d.m. D\ phrases of the type "D\ is isolated" are used throughout the chapter. 



3.1 Possible configurations 

The aim is now to determine all the possible 2s-cycle configuration (using (1)- 
presentations) that can arise in a matrix. Some new terminology and lemmas are 
introduced first. 

Definition 3.11. The weights vector of a 2s-cycle configuration is defined as the 
vector containing all the tij 's in such configuration. 
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For example the cycle configurations presented previously: 



V = 2 2 



V = 



2 2 
2 



V = 



2 2 
1 
1 




1 
1 



have weights vectors 



[2,2], [2,2,2], [2,2,1,1,1,1]. 



Theorem 3.12. The weights vectors that can make a 2s-cycle configuration are 
all the possible combinations of numbers, that sum to 2s and that do not contain 
exactly two odd values. 

Proof. The sum of the elements in a weights vector must be 2s for the definition. 
Hence all the possible vectors that sum to 2s are candidate to be weights vectors of 
a 2s-cycle configuration. The vectors that have two odd weights can be discharge 
since in such cases it is impossible to have even row weight and column weight for 
each d.r. and d.c. In fact, if both sub-matrices with odd weights are in the same 
d.r. the two d.c.'s where they lie have odd weight, and vice-versa. □ 

The following theorem specifies which of these weights vectors can form a 
cycle configuration of chosen dimension. 

Theorem 3.13. An (r, c) 2s-cycle configuration can be obtained only from weights 
vectors that have the following characteristics: 

• they contain at least r + c — 1 elements, 

• they contain no more than min(cr, 2s) elements, 

• the value of the biggest element is no more than 2s — 2(c — 1) 

Proof. Suppose r < c, for remark 13.71 the case r > c can be reduce to this by 
transposing the configuration. Any d.r. must have at least one element with tij ^ 
otherwise it would be a (r — 1, c) cycle configuration, the same is true for the d.c.'s. 

The value ti t % must not be for the definition of (l)-presentation, then for 
Lemma [3.101 there must be at least another d.m. with / in the same d.r. or 
d.c. Suppose, without loss of generality, that tip ^ 0. Applying Lemma [3. 101 to 
the sub-matrix [Ti i, Tip] implies that there must be another d.m. with tij ^ in 
the first d.r. or in one of the first two d.c.'s, otherwise [Ti i , Tip] would be isolated. 
Assuming ti^ ^ another sub-matrix, [Ti i, Tip, ?i 3], is obtained. Repeating 
the process until all t\ x 7^ 0, r — 1 d.r.'s that must have at least one tij 7^ 
element are left, hence a total of at least c + r — 1 weights are needed. Note that 
the same result will be obtained if instead of "filling" the d.c.'s first the d.r.'s are 
filled, or any combination. It is obvious that any weights vector cannot have more 
that 2s values. Moreover, the sums of the elements in a weights vector is 2s but 
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there cannot be more elements that d.m. hence the max number of elements in a 
weights vector is min(2s, rc). Using the (l)-presentation then ti t i > tij for every 
i, j, the remaining d.c. must not be zero and have at least row weight 2. Since the 
sum of all row weights must be 2s, the max value of t% \ is 2s — 2(c — 1). □ 



Using theorems 13.121 and 13.131 it is possible to determine the set of possible 
weights vectors that can be used to form a (r, s) 2s-cycle configuration. Obtaining 
the configurations from the weights vectors is a mater of placing the weights in 
such a way that they satisfy Lemma [3^91 and Lemma 13.101 

Next, an example is presented. The aim is to find which configurations of type 
(4, 5) can have 10-cycles. The process starts by looking for all possible vectors of 
number that sum to 10. 

[10], 



[8,2], jp^, 

)m(,j^?T, [7,i,i,i], 

[6, 4] , pS^, [6, 2, 2] , Jg^<lf , [6, 1, 1, 1, 1] , 

J^lf , J^^f, [5, 3, 1, 1] , J^MCEf, [5, 2, 1, 1, 1] , [5, 1, 1, 1, 1, 1] , 
[4, 4, 2] , Ta^M^f, T^i^f, J^MClf , [4, 3,1,1,1] , [4, 2, 2, 2] , J4^<lf, 

[4,2,1,1,1,1] ,[4,1,1,1,1,1,1] , " ' 
[3, 3, 3, 1] , J^<5f , [3, 3, 2, 1, 1] , [3, 3, 1, 1, 1, 1] , 1^2?%^, t 3 ' 2 ' 2 ' ^ ' 

[3, 2, 1,1, 1,1,1], [3, 1,1, 1,1, 1,1,1], 
[2, 2, 2, 2, 2] , jgj^^TXlT, [2, 2, 2, 1, 1, 1, 1] , [2, 2, 1, 1, 1, 1, 1, 1] , 

[2, 1,1, 1,1, 1,1, 1,1], [1,1, 1,1, 1,1, 1,1, 1,1]. 

The crossed out vectors are not to be considered since they contain exactly two odd 
numbers. Applying Theorem 13. 131 it is possible to eliminate weights vectors that 
do not have length between 8(= 5 + 4 — 1) and 10(= min(10, 20)) and that do not 
have max value less or equal to 2(= 10 — 2(5 — 1)). The remaining candidates are 

[2, 2, 1, 1, 1, 1, 1, 1] , [2, 1, 1, 1, 1, 1, 1, 1, 1] , [1, 1, 1, 1, 1, 1, 1, 1, 1, 1] . 

It is now necessary to place the values of the weights vectors inside the (4, 5) d.m.. 
It is now proved how the only possible configurations are : 



1. 



2 110 

10 10 

10 1 

1 1 



11110 
1 1 
110 
1 1 



This can be proved with some easy considerations. There are twenty Tjj with 

1 < i < 4, 1 < j < 5 . It must be £\ . t^j = 10, and t^i > tij, Vi, Vj s.t. i < 
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4, j < 4 . There are five columns and for any of this column the weight must be 
even and not zero, (if there is a zero column then it falls in a smaller configuration), 
the only possibility to have a total sum of ten and column weights even is that all 
columns have weight 2. Moreover, there are four rows and for any of these rows 
the weight must be even and not zero, (if there is a zero column then it falls in 
a smaller configuration), the only possibility for the total sum to be ten and row 
weights even is to have one row with weight 4 and the remaining with weight 2. 

• Considering the weights vector [2,2,1,1,1,1,1,1]. 

An element ti j = 2 cannot be on a d.r. of weight 2 because both the d.c.and 
d.r. would be completed but this cannot be otherwise it would be isolated ( 
Lemma I3.10I ). So, both the elements t^j = 2 must be in the same row, the 
only row with weight 4. In this case such a row has weight 4 and cannot 
have other non zero elements in it, but also the two columns are completed, 
since they have weight 2. This situation does not satisfy Lemma l3.10l hence 
this weights vector does not lead to any cycle configuration of this size. 

• Considering the weights vector [2,1,1,1,1,1,1,1,1]. 

Element ti t % = 2 since it is the max. Lemma 13-101 implies that there must 
be at least another non zero element in row one, but since only weight of 
value one are present and the row weight must be even then there must be 
two one elements in row one. Supposed that ti^ = ii,3 = 1. Considering 
now Lemma [3791 applied to the remaining rows and columns configuration 1, 
or a column/row permutation of it, is found. 

• Considering the weights vector [1, 1,1,1,1,1,1,1,1,1]. 

It can be supposed that the first row has weight 4 hence it must have four 
elements tij = 1 that can be ti i = ti 2 = £1,3 = ii,4 = 1. Applying 
Lemma 13.91 to all remaining rows and columns configuration 2, or a col- 
umn/row permutation of it, is obtained. 

And this prove the claim. 



3.2 4-cycle configurations for generic H matrices 

In this section all the possible configurations that can give cycles of length 4 in the 
case of a generic matrix are found and listed. 

Theorem 3.14. Let B 6 M m m,/3,7- The only possible 4-cycle configuration, in 
(l)-presentation, are as follows^; 

1. (1,1), 

|4|, 

2- (1,2), 

|22|, 
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3. (2,2), 



1 1 
1 1 



Proof. Following from Theorem 13.121 the only possible weights vectors for a 4- 
cycles are: 

[4], [2,2], [1,1,1,1]. 

From them it is straightforward to find the configurations listed in the statement. 

□ 



3.3 6-cycle configurations for generic H matrices 

In this section all the possible configurations that can give cycles of length 6 in the 
general case are found and listed. 

Theorem 3.15. Let B € M.m,a,3a- The only possible 6-cycle configurations, in 
(l)-presentation, are as follows 2 : 



1. (1,1), 

2. (1,2), 

3. (1,3), 
4- (2,2), 

5. (2,3), 

6. (3,3), 



|6[, 
| 4 2 |, 
2 2 2 I 



4.1 



2 2 
2 



4.2, 



3 1 
1 1 



2 11 
Oil 



110 
1 1 
Oil 



For brevity any configuration that is the transpose of another is omitted. 
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Proof. Lemma [3761 proves that these are all the dimensions that 6-cycle configura- 
tions can have. It is now necessary to prove that the listed configurations are the all 
and only possible configuration that can have 6-cycles. To do so it is necessary to 
prove that for each dimension all the possible configurations associated it are listed 
in the statement. 

Following from Theorem 13. 12l the only possible weights vectors for a 6-cycles 



The weights vectors are used to simplify the process of determining the valid 
configurations. The study of each case is presented next. 

Configuration (1,1). A type (1, 1) configuration can be generated only by weights 
vector [6] and it corresponds to caseQ] 



Configuration (1,2). A type (1,2) configuration can be generated only by weights 
vector [4, 2] and it corresponds to case [2] 



Configuration (1,3). A type (1,3) configuration can be generated only by weights 
vector [2, 2, 2] and it corresponds to case [3] 

Configuration (2,2). Applying Theorem 13. 131 configurations of type (2,2) must 
have one of the following weights vectors: 



Applying Lemma [3~9l it is straightforward that these two weights vectors result in 
configurations [3] 1 and[3]2. 



Configuration (2,3). Applying Theorem 13. 131 configurations of type (2,3) must 
have one of the following weights vectors: 



• [1, 1, 1, 1, 1, 1] is not a possible choice because it would require to have an 
odd (3) row weight since there are only two d.r.'s 

• [2, 1, 1, 1, 1] results in configuration [5] This can be proved with some short 
considerations. The value t\ \ = 2 since it is the max value and £2,1 = 
because otherwise column weight of c-1 would be odd. It follows that all the 
other must be 1 and this proves the claim. 



are: 




[2, 2, 2], [3, 1,1,1]. 



[2, 1,1, 1,1], [1,1, 1,1, 1,1]. 
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Configuration (3,3). Applying Theorem 13. 131 configurations of type (3,3) must 
have one of the following weights vectors: 



[2, 1,1, 1,1], [1,1, 1,1, 1,1] 



• [2, 1, 1, 1, 1] is not a possible choice. In fact, there are three d.r.'s and for 
any of this d.r.'s the row weight must be even and not zero (if there is a zero 
row then it fall in a smaller configuration). The sum of all d.r.'s must be six, 
hence all d.r's must have row weight 2. The same is true for the d.c.'s. If any 
tij = 2 then that element is isolated (Lemma [3. 101 ) but this is not allowed 
hence this weights vector cannot generate valid configurations. 

• [1,1,1,1,1,1] results in configuration [6] To prove this it is sufficient to con- 
sider that as discussed in the previous configuration each d.r. and d.c. mush 
have weight two. Any possible way to put two 1-elements in each d.r. and 
d.c, satisfying Lemma 13.101 results in configuration [6] or in a column/row 
permutation of it. 

Thanks to Remark 13.71 it is not necessary to prove the transposed configurations 



3.4 8-cycle configurations for generic H matrix 

In this section all the possible configurations that can give cycles of length 8 in the 
case of general matrices are found and listed. 

Theorem 3.16. Let B € -Mm,^^- The only possible 8-cycle configurations, in 



□ 




1. (1,1), 



8 



2. (1,2), 



2.1 6 2 



2.2 4 4 



3. (1,3), 



4 2 2 



4. (1,4), 



2 2 2 2 



3 For brevity any configuration that is the transpose of another is omitted. 
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5. (2,2), 



5 


1 


, 5.2 


4 


2 


, 5.3 


4 


2 


1 












1 


2 


2 


3 


1 


, 5.5 


3 


1 


, 5.6 


2 


2 




1 


1 








3 


3 


2 


2 



6. (2,3), 



6.1 



6.4 



5.1 



5.4 



4 1 1 

1 1 

2 2 2 

2 



6.2 



6.5 



3 2 1 
1 1 



1 1 
1 1 



6.3 



6.6 



3 10 

1 1 2 

2 2 
2 2 



7- (2,4), 



7.1 



2 2 11 
11 



7.2 



2 110 
112 



7.3 



1111 
1111 



(3,3), 





CO 


1 







2 


1 


1 




2 


1 


1 




2 








1 


1 





1 


, 8.2 


2 








, 8.3 


1 


1 





, 8.4 


1 


1 










1 


1 







1 


1 




1 





1 




1 


1 


2 



9. (3,4), 



2 110 
10 1 
11 



9.2, 



1111 
110 
11 



10. (4,4), 



1 


1 








1 





1 








1 





1 








1 


1 



Proof. The proof of this theorem is long and similar to the proof of theorems 13.141 
and 13.151 with the addition of some logical reasoning. To improve readability 
it is here omitted, an interested reader can find the full and detailed process in 
Appendix. 

□ 
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The previous three theorems list all the configurations that can contain cycles 
of length less than ten. In particular such configurations identify how many cycle 
points lie in each sub-matrix (d.m.). To improve such result it is necessary to endow 
the d.m. with some structure that allows to predict the positions of such point. The 
next section considers the case where the d.m.'s are circulant matrices. 

4 The quasi-cyclic case 

This section restricts the discussion to matrices which can be used as parity-check 
matrices for LDPC quasi-cyclic codes, and gives some generic definitions and 
lemmas. 

A quasi-cyclic code of index t is a linear block code C in which a cyclic shift 
of any codeword in C by t positions is also a codeword. The generator matrix G 
for these codes is a matrix where every row is a t circular shift of the previous 
row. It can be shown that every generator matrix for a quasi-cyclic code can be de- 
composed in circulant matrices. This property makes the encoding suitable for low 
cost hardware implementation. Major results on a Array codes, a class of quasi- 
cyclic codes, have been presented by Fan 11311 and later by Fossorier ll28l . Array 
codes are composed by J rows of L circulant matrices. The circulant matrices used 
was limited to have weight- 1. The condition given by Fan [31] (later re-proposed 
in Theorem 2.1 in [28]) gives the conditions for the existence of cycles in such 
class and is rewritten here using the notations presented previously. We remind the 
reader that with our notations e(j>) is any of the exponents of the polynomial p and 
s(p) is the separation of the two monomials. 

Theorem 4.1 (Fan-2000). A necessary and sufficient condition for the matrix H to 
have a Is-cycle is: 



k=0 

were e(pi and e(j>2 k) are me exponents of the circulant matrices that contain 
the two cycle-points that lie in the same cycle-column k. 

Remark 4.2. The theorem was given in the case of a particular case of quasi- 
cyclic codes with only weight 1 circulants but the same theorem holds in the case 
when weight-2 circulant matrices C are present. The presence of weight-2 circu- 
lant matrices allows cycle columns to lie in the same circulant. In such case the 
difference e(pik) ~ € (P2,k) JJ equivalent to the s(pk) of such circulant. Previous 
work considering weight-2 circulant matrices was presented by Smarandache and 
Vontobel et al. [32] but focuses on the problem of minimum distance and not girth 
properties. 

However, the construction by Fossorier cannot provide codes whose Tanner 
graph has a girth higher than 12, as shown by Fossorier himself in the same paper 
and known from Fan. To overtake such limitations and to complete the study of 



s-1 




22 



the cycles in quasi-cyclic codes a wider class of quasi-cyclic LDPC codes are 
considered in this thesis. 



Definition 4.3. Let B G M. m a,8 7 an d let {A{j} form its standard decomposition. 
The matrix B is said to be in C ma p~ if any d.s. A^j can only be either a weight-2 
circulant matrix, a weight-1 circulant matrix or an m x m zero matrix. 

It is well-known that any quasi-cyclic code has a parity-check matrix H made of 
circulant sub-matrices. However, if an LDPC code with good girth characteristics 
is desired, it is necessary to avoid sub-matrices that are weights circulant matrix, 
with t > 3, since they contain internal 6-cycles [33 ] and codes with higher girth 
are wanted. Matrices in C TO|Q)( g !7 are the only interesting parity-check matrices for 
good quasi-cyclic LDPC codes. This family of H matrices is the focus of the study 
presented here. 

In the case of quasi-cyclic LDPC codes a cycle configuration c in C m ai/ g i7 will 
be described as: 



(7) 



C-2 J-l C-l 
O J-l J-l 



where C — 2 is a weight-2 circulant matrix containing two points of the cycle, O is 
a zero matrix, J — 1 is a weight-1 matrix containing one point of the cycle, C — 1 
is a weight-2 circulant matrix containing one point, and so on. Clearly, notations 
of the type O — i are not used because it is obvious that in a zero matrix there must 
be i = 0. In a situation like (0 it will be written that c "contains" a matrix (d.s.) 
C — 2, three matrices J — 1 and a matrix C — 1. When no confusion can arise, even 
looser notation will be adopted. For example, in ([7J) it may be said that c contains 
C — 2, that c contains C — 1 and that c contains J — 1. Similarly, phrases of the 
type "c contains d.m.'s, d.s.'s,d.r.'s, d.c.'s" will be used with the obvious meaning. 
For example, in (0) c contains the following d.m.'s 



C-2 J-l C-l , \ O J-l J-l 



C-2 J-l 
O J-l 



It is sometimes convenient to gather together cycle configurations possessing a 
given number of cycle points. To be more precise, the notation A — i means that 
A can be any allowable matrix containing i points. In other words, the notation 

I C-2 A-2 I , 

is equivalent to the union of 

\ C-2 C-2 I and \ C-2 J-2 \ . 

Let c be a cycle configuration, the following easy lemmas can be used to show 
that c cannot exist in C m ^ a ^ n . The expression "to discard configuration c" is syn- 
onymous with "to show that configuration c cannot exist in C m ^ a ^ a ". 
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Lemma 4.4. Suppose c contains J — i, with i > 2. Then 

1. J — i can contain neither a cycle column nor a cycle row, 

2. in c, any d.r. and any d.c. containing J — i (i > 1) must contain another 
non-zero matrix, 

3. c can be discarded if it contains only one d.r. or only one d.c. . 

Proof. Point one is obvious, because a J is a shift of an identity matrix hence there 
is only one non zero entry for each row and column. Point two follows from point 
one and the fact that the cycle point contained in the J must be part of a linked set. 
Part three follows from part two. □ 

Lemma 4.5. Ifc correspond to a 2s-cycle with 2s < 8, it cannot contain a weight-1 
circulant matrix J — i with i > 3. 

Proof. Otherwise any cycle point in J will need another point in the same row and 
one in the same column, to form a linked set. For lemma l4~4l this points cannot be 
in the J but this implies the existence of at least nine points in the linked set, but 
this defies the definitions of 2s-cycle with 2s < 8 □ 

The following two lemmas are obvious and are reported here only for clarity. 

Lemma 4.6. Suppose c has a (1,2) d.m. composed of two weight-1 circulant 
matrices 

| J — i J — 1\, i > 1. Then in that d.m. there cannot be more than one cycle row. 
Similarly for a (2, 1) d.m. in c. 

Lemma 4.7. Suppose that c contains a matrix J — i, i > 1. If C — i is substituted 
to J — i, another possible cycle configuration is obtained. 

Passing from a cycle configurations obtained with the procedure outlined in the 
previous section to the respective configurations in the quasi-cyclic case require 
the application of lemma l4~4l lemma l4~6l lemma 14771 and lemma |4~51 to the original 
configuration. For example, configuration: 



2 110 

10 10 

10 1 

1 1 



results in a cycle configuration for quasi-cyclic codes 



C-2 A-l A-l 

A-l 

A-l 




A-l A-l 





1 

A-l 
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In fact cycle configuration: 



J-2A-1A-1 
A - 1 1 

A-l A- 1 ' 

A-1A-1 

can be discard because in J — 2 there is no cycle column (lemma l4~4l . 

In the following subsections the theorems presented previously (13.14ll3~T5ll3.16l) 
are specialized for the quasi-cyclic case. All the cycle configurations that can ap- 
pear in such case are listed. For every circulant in a configuration the possible 
weights that such a circulant can have is discussed. 

4.1 4-cycle configurations for quasi-cyclic H matrices 

Theorem 4.8. Let be M G C m ,a,/3,7- The configurations in M that may contain 
cycles of length 4, are the following: 

1. 

|C-4|, 

2. 

\C-2 C-2|, 

3. 

A-l A-l 
A - 1 A - 1 ' 

Proof. All the configurations, for the generic case, that appeared in theorem f3. 141 
For each configuration it is proved how only configurations listed in the theorem 
are valid in the case of a quasi-cyclic LDPC code. 

Configuration [J gives 

|C — 4] . 

Since cycle configuration | J — 4| may be discard (Lemma 1431) . 
Configuration |2] gives 

| C-2 (7-2 | . 

Cycle configurations \C — 2 J — 2| and | J — 2 J — 2[ may be discard because 
in J — 2 there is no cycle column (Lemma l4~4l -l). 



4 For brevity any configuration that is the transpose of another is omitted. 



25 



Configuration [3] gives 



A-l A-l 
A-l A-l 



Note how, thanks to Lemma [4771 it is necessary only to show that the following 
configuration is possible: 



and this is obvious. For Remark 13771 it is not necessary to study the cycle configu- 
rations that are transposed of the one considered. Hence it has been proved that the 



To better explain the meaning of A and what it implies the (non-equivalent) 
cycle configurations present in the statement of Theorem 14. 8 l ease l3l , are reported: 



J - 


1 


J - 


1 




c- 


1 


J - 


1 




C-l 


C-l 


J - 


1 


J - 


1 


1 


J - 


1 


J - 


1 




J -I 


J -I 


c- 


1 


J - 


1 




c - 


1 


c - 


1 




C-l 


C-l 


c - 


1 


J - 


1 


1 


c - 


1 


J - 


1 


? 


C-l 


C-l 



4.2 6-cycle configurations for quasi-cyclic H matrices 

Theorem 4.9. Let be M E C m]Qi/ g A . The configurations in M that may contain a 
cycle of length 6, are the following^ 



J -I J -I 
J -I J -I 



listed configurations are the only valid ones. 



□ 



1 



C-6 



2. 



C-4 



C-2 



3. 



C-2 C 



2 C-2 



4. 



C-2 
O 



A-2 
C-2 



5. 



C-2, J-l 
J -I J-l 



5 For brevity any cycle configuration that is the transpose of another is omitted. 
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6. 

C-2 A-l A-l 

O A-l A-l ' 

7. 

A-l A-l O 

A-l O A-l . 

O A-l A-l 

Proof. All the configurations that appear in Theorem 13.151 are considered and it 
is proved how there exist no other cycle configurations from the one listed in the 
statement. The process also proves that no unnecessary cycle configurations are 
listed. 

Configuration Q] gives 

|C-6| . 

Since cycle configuration j J — 6| can be discard (Lemma |43T ). 
Configuration |2] gives 

| C-4 C-2 | . 

Other cycle configurations |C — 4 J — 2|, | J — 4 J — 2| and | J — 4 C — 2| may be 
discard, because in J — 2 and in J — 4 there is no cycle column. 

Configuration |3] gives 

| C-2 C-2 C-2 | . 

Other cycle configurations may be discard because they contain J — 2 and in 
J — 2 there is no cycle column (Lemma l4~4l -1). 

Configurations] 



1. Configuration |Hl gives 



C-2 C-2 
O C-2 



C-2 J-2 
O C-2 



Other cycle configurations may be discard because they contain J — 2 as the 
only non-zero matrix in a d.r. or d.c. (Lemma l4~4l -2). 



2. Configuration @1 2 gives 



C 
A 



3 A-l 
1 A-l 
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In fact any cycle configuration of kind 



J-3 A - 1 
A-l A-l 



can be discarded (lemma 14751 ) . On the other hand, cycle configuration 

C-3 J-l 
J-l J-l 



is obviously acceptable and so lemma 14771 can be applied. 
Configuration |5] gives 



C-2 A-l A-l 
O A-l A-l 

In fact , the following cycle configuration is obviously possible 



C-2 J-l J-l 
O J-l J-l 



Hence lemma l477l can be apply. 

To discard all cycle configurations of type 



J-2 A 
O A 



1 A-l 
1 A-l 



it is enough to apply Lemma 14741 -2 to d.c. 1. 
Configuration gives 

A-l A-l O 
A-l O A-l 
O A-l A-l 

Clearly, the following cycle configuration is possible 



J-l J-l O 
J-l O J-l 
O J-l J-l 



hence Lemma[4/7jcan be applied. 

For remark 13.71 it is not necessary to study the cycle configurations that are 
transposed of the one considered. Hence it has been proved that the listed configu- 
rations are the only valid ones. □ 
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4.3 8-cycle configurations for quasi-cyclic H matrices 



Theorem 4.10. Let be M € C m)OL fi^. The configurations in M that may contain 
a cycles of length 8, are the following^ 



1. 



3. 



|<7-8|, 
|C-6 C-2|, 
|C-4 <7-4|, 



1(7-4 C-2 C - 21 . 



1(7-2 C-2 (7-2 (7-21 



7. 



(7-5 A-l 

A-l A-l 

(7-4 A-2 

C-2 



(7-4 (7-2 
(7-2 



10. 



11. 



12. 



(7-3 (7-3 

A-l A-l 

(7-3 A-l 

A-l (7-3 

A-2 A-2 

A-2 A-2 



(7-4 A-l A-l 
O A-l A-l 



For brevity any cycle configuration that is the transpose of another is omitted. 
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13. 



14. 



15. 



16. 



17. 



18. 



19. 



20. 



C-3 C-2 A-l 

A-l O A-l 

C-3 O A-l 

A-l C-2 A-l 

A-2 C-2 C-2 

C-2 O O 

A-2 A-l A-l 

A-2 A-l A-l 

C-2 A-2 O 

O A-2 C-2 



C-2 C-2 A-l A-l 

O O A-l A-l 

C-2 O A-l A-l 

O C-2 A-l A-l 

A-l A-l A-l A-l 

A-l A-l A-l A-l 



21. 



22. 



23. 



C-3 A-l O 

A-l O A-l 

O A-l A-l 

A-2 A-l A-l 
C-2 O O 

O A-l A-l 

A-2 A-l A-l 

A-l A-l O 

A-l O A-l 
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C-2 O O 

A-l A-l O 

A-l A-l C-2 

25. 



26. 



27. 



c - 


2 


A - 


1 


A- 


1 












A - 


1 







A- 


1 







O 




A- 


1 


A- 


1 


A- 


1 


A - 


1 


A - 


1 


A- 


1 


A- 


1 


A - 


1 






















A - 


1 


A- 


1 


A- 


1 


A - 


1 












A- 


1 







A - 


1 


o 









A - 


1 







A- 


1 












A - 


1 


A- 


1 



Proof. Once more for brevity the full proof is omitted here and can be found in 
Appendix. □ 



5 Relations between polynomials and the existence of cy- 
cles 

This subsection gives an interpretation to the cycle configurations presented pre- 
viously in terms of the polynomials associated to the circulant matrices involved. 
First a general statement on the girth of the Tanner graph associated to a binary 
weight-2 circulant matrix is given. The following proposition has first been pre- 
sented in ll34l and |[35t . Here only the main result is reported. 

The notation introduced in subsection l2.3l is presented again here for clarity: 

• e{p) : any of the exponents (a, b) of the polynomial p(x) = x a + x b , 

• s(p) : the separation (min(6 — a, a + m — b)) of the polynomial p. 

Proposition 5.1. Let m > 3. Let M = (Mjj) be a circulant binary mxm matrix, 
generated by a weight-2 polynomial p. Let s(p) be the separation ofp and g be the 
girth of the Tanner graph of M. Then 



gcd(m, s(p)) 
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Proof. Let G be the Tanner graph of M. Then each check node is connected 
exactly to two bit nodes, and vice-versa. We denote by c\, . . . ,c m and b±,..., b rn , 
respectively, the check nodes and the bit nodes of G. Without loss of generality, 
we may suppose that = 1. By circularity, Mjj = 1 for 1 < j < m, which 
implies that cj is connected to bj, for 1 < j < to. 

By definition of separation, we have either s = b — a or s = m — (b — a). 
By circularity we may assume s = b — a, so that Mi iS +i = 1 and b — a < to/2. 
But then we have Mj :S+ j = 1, for 1 < j < m — s, and M ?iS+ j_ m = 1, for 
to — s + 1<j'<to. Therefore, check node Cj is also connected to: either bit node 
b s+ j, which happens when 1 < j < to — s, or bit node b s+ j- m , which happens 
when to — s + 1<j'<to. No other non-zero entry is present in M and so no 
other connection exists among nodes in G. 

For any i, we want to find the length gi of the minimum length cycle containing 
Q. The girth of G will be g = mini<j< m g^. By symmetry of G, gi does not depend 
on i, in particular g = g\. 

We now determine g\, constructing a path as follows: 

• we start from c±. From c\, we may go either to b\ or to b s +\. We choose to 
go to b\. Clearly, the cycle will be closing when we will find ourself in b s+ i, 
as the next step will be c\ . From now on, the path allows no more choices. 
We will use arrows to shorten our notation. 

• ci — > b\, b\ — > c m _ s+ i. If to — s + 1 = 1, the cycle will be closed and 
gi = 2: this is impossible since s < to. 

• c m - s+ i — > b rn - s+ i, b m - s+ i — > c m _2s+i (to — 2s + 1 > 1). 
If to — 2s + 1 = 1, the cycle is closed and g\ = 4. 

• We perform steps of type 

c m~(l-l)s+l ~^ b m -(l-i) s+ i, & m _(7_i) s+ i — > C m _; s+ i , 

until either there is an / > 1 s.t. to — Is + 1 = 1 or there is no such I. We 
analyze the two cases separately. 

There is an I s.t. to = Is. This is equivalent to s\m. This is also equivalent to 
<7i = 21 because we have formed a length-2Z cycle and we have encountered 
no smaller cycles. 

Or there is no / s.t. m — Is + 1 = 1. In this case let I = \m/s\. Then g > 21, 
i.e. g > 2m/ s. The next two steps will be 

c m-ls+l ~^ b m _i s+1 , b m _i s+1 — ► C 2m _(7 +1 ) s+1 , 

i.e. we have to "wrap on the graph". 

• We perform steps of type 

c 2m-(7-l)s+l —> &2m-((-l)s+li &2m-(7-l)s+l ~* c 2m-ls+l , 
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until either there is an I > \m/s\ s.t. 2m — Is + 1 = 1 or there is no such I. 
We analyze the two cases separately. 

There is an I > \m/s\ s.t. 2m = Is (but no I > 1 is s.t. m = Is). Obviously 
this is equivalent to s\2m and not s\m. 

This is also equivalent to g\ = 21 because we have formed a length-2/ cycle 
and we have encountered no smaller cycles. 

Or there is no I s.t. 2m — Is = (but no I > 1 is s.t. m = Is). This is equiv- 
alent to s not dividing 2m (and s not dividing m). This is also equivalent to 
g > \2m/s\. In this case let T= \2m/s\. The next two steps will be 

c 2m-(s+l ~~ * ^2m-ls+li ^2m-ls+l ~~ ^ c 3m-(r+l)s+l ' 

i.e. we have to "wrap on the graph" again. 

• In the general case, after 21 steps we reach c gm _; s+ i, where s does not 
divide any zm for 1 < z < q— 1, by a generalization of the above arguments. 
At some stage we must meet again c%, so that c gTO _/ s+ i = c\. In other words, 
the girth is 21, with I = qm/s, if q is s.t. 

s\qm, s J(zm, 1 < z < q — 1 . 

This means that qm is the smallest multiple of m that is also a multiple of s, 
i.e. qm is the minimum common multiple of m and s. 
Let [, ] denote minimum common multiple. For any two integers a and b, we 
have [a, b]/a = bj gcd(a, b), so that 

qm [m, s] m 
s s gcd(m,s) 

□ 

Remark 5.2 (Minimum dimension of the circulants). Theorem \5. 1 1 indirectly gives 
a minimum dimension that the circulant matrices must have to allow high girth. In 
particular ifH matrices with girth 10 (or higher) are desired the circulant matrices 
must be at least [5 x 5]. To prove this is sufficient to consider that to have g > 10 

it must be — ,, m , > 5, hence m > 5. 

gcd(m,s(p)) — — 

A plot of the girth that can be obtained given a certain circulant dimension, 
changing the value of the separation, is presented in Figure [3] It is evident that the 
maximum girt (g = 2m) is also achievable, e.g. by choosing the separation to be 
one. Of more interest is to consider that which values of m and s must be avoided 
to guarantee a good girth. 

The following theorems present which conditions, on separations and expo- 
nents, must be satisfied for a certain cycle, in a cycle configurations, to exist. Once 
this set of conditions is known matrices with high girth can be generated by choos- 
ing the polynomials of the circulants in such a way that none of the conditions are 
satisfied. 
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3/-£ 10 11 12 IS- 14 15 16 17 18 19 20 

Circulant dimension m 



Figure 3: Girth values that can be obtained with a single circulant given its dimen- 
sion 

5.1 Conditions for the 4-cycles 

Theorem 5.3. Let be M E C m)aj/ g i7 . The configurations in M that may contain a 
cycles of length 4 and associated conditions, are the following 

1. 

\C - 4| , s(p) = m/2 

2. 

\C l -2 C 2 -2\, s(p 1 ) = s(p 2 ) 

3. 

e (pl)_ e (p2)_ e (p3) +e (p4 )=0 

Proof. It is of interest to note that the conditions associated with every configu- 
ration can be proved applying theorem 14- 1 1 that is the generalization of Fossorier 
theorem for weight-2 circulants. 

It has been chosen to follow a different methodology that slightly reduces the 
notation and allow to use a graphical representation that easier to understand. The 
use of separations instead of the difference of exponents makes the conditions 



A 1 -! A 2 -l 
A 3 -l A 4 -l 
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clearer and more evident; it also reduces the number of variables involved mak- 
ing them easier to check. The proof is divided in shorter lemmas each considering 
a case of the main theorem. An 4-cycle is formed by two cycle columns called y, z 
and two cycle rows called x, t. In every 4-cycle there are 4 cycle points, they take 
the name of the column and row where they lie in (x, y), (x, z), (t, y), (t, z). The 
notation is used to compute the conditions. To allow the reader to better understand 
the meaning of the formulae for each lemma a graphical representation of the cycle 
is presented. For example Figure @]presents an example of 4-cycle on two weight-2 
circulants. The two rectangular blocks represent the two circulants and the dashed 
diagonal lines the position of the ones in the matrices. The cycle columns and cycle 
rows are marked by the dash-doted lines and the cycle is outlined with a wider line. 

Lemma 5.4. There is a 4-cycle in case\l\if and only if 

2\m and s(p)=m/2. 

Proof. There is a 4-cycle if and only if g = 4, because smaller girths are not 
possible. Applying Proposition 15.11 to the case 5 = 4 and M = C, there is a 
4-cycle if and only if 

m 

4 = 5 = 2—- 

gcd(m, s) 

i.e. m = 2gcd(m, s). In particular, m is even and m/2\s, but s < m/2, so that 
s = m/2. On the other hand, if m is even and s = m/2 then gcd(m, s) = s = 
m/2. □ 

Lemma 5.5. There is a 4-cycle in case\2\if and only if 

sip 1 ) = s(p 2 ) 

Proof. It can be assumed that there is a 4-cycle if and only if column y lies in C 1 
and column z lies in C 2 (Figure [4]). 



Figure 4: Example of 4-cycle on two weight-2 circulants. 
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Applying Prop. 12.91 -2 to column y yields 

(8) x-t = ±s 1 . 
Applying Proposition 12.91 -2 to cycle column z yields 

(9) x-t = ±s 2 . 

From © and ©, s 1 = ±s 2 is obtained and hence s 1 = s 2 (Lemma [2.10l ). □ 

Lemma 5.6. There is a 4-cycle in case\3\if and only if 

e(p l )-e( P 2 )-e(p 3 )+e(p 4 ) = 

Proof. It can be assumed that there is a 4-cycle if and only if, simultaneously, cycle 
point (x, y) lies in A 1 , cycle point (x, z) lies in A 2 , cycle point (t, y) lies in A 3 
and cycle point (t, z) lies in A 4 (Figure [5]). 
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Figure 5: Example of 4-cycle on four weight- 1 circulants. 
Since cycle point (x, y) lies in A , using Proposition 12.91 1: 
(10) y = x + eip 1 ) . 
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Since cycle point (x, z) lies in A 2 , using Proposition 12.91 -1: 



(11) z = x + e(p 2 ). 
Since cycle point (t, y) lies in A 3 , using Lemma l3~4l 

(12) y = t + e(p 3 ). 
Since cycle point (t, z) lies in A 4 , using Lemma l3~4l 

(13) z = t + e(p 4 ) . 

The desired result is obtained from ([TOl), £EU>, GSJ) and (fT3l) . 



□ 



It has hence been proved that the conditions listed on the statement considers 
all the possible 4-cycles that can exist on the studied quasi-cyclic matrices. □ 

For example consider the polynomial p(x) = 1 + x 3 with m = 6 for such 
polynomial s(p) = 3 hence s(p) = m/2 and for condition 1 a 4-cycle exist. 



1 
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1 1 
5.2 Conditions for the 6-cycles 

Theorem 5.7. Let be M G C m a p ~. The configurations in M that may contain a 
cycles of length 6, are the following 



1. 



2. 



3. 



IC-61 



|C 1 -4 C 2 -2|, 



s(p) = m/3 



s(p 2 ) = ±2*(p 1 ) 



|C 1 -2 C 2 -2 C 3 -2|, sCp 1 ) ±s(p 2 ) ±s(p 3 ) = 0, 



37 



4. 



C l - 2 A 2 - 2 
O C 3 -2 



e(p 2 )-e(p 2 ) = ±s(p 1 )±s(p 3 ), 



C 1 - 3 A 2 - 1 
A 3 -l A 4 -l 



e (pl)_ £ (p2)_ € (p3 )+€(p 4 )s±a(pl)j 



C 1 -2 A 2 -l A 3 -l 
O A 4 -l A 5 -l 

e(p 2 )-e(p 3 )-e(p 4 ) + e(p 5 ) = ±s(p 1 ), 



A 1 - 1 A 2 - 1 O 
A 3 - 1 O A 4 - 1 
O A 5 - 1 A 6 - 1 

e(p l) _ e(p 2 ) _ e (p3) + e (p4) + _ e( /) 



Proof. As for the previous theorem the proof is once more divided in smaller 
lemmas each considering a particular case. An 6-cycle is formed by three cy- 
cle columns called y,z,v and three cycle rows called x,t,w. In every 6-cycle 
there are 6 cycle points, they take the name of the column and row where they 
lie: (x,y), (x,z), (t,y), (t,v), (w,v), (w,z). The notation is used to compute the 
conditions. 

Lemma 5.8. There is a 6-cycle and there are no A-cycles in case\l\if and only if 

2>\m and and s(p)=m/3. 

Proof. There is a 6-cycle and no 4-cycle if and only if g = 6. Applying Prop. 15.11 
to the case g = 6 and M = C, this is equivalent to 



gcd(m, s) 

i.e. m = 3gcd(m,s). In particular, 3|m and m/3\s, but s < m/2, so that s = 
m/3. On the other hand, if m is divisible by 3 and s = m/3 then gcd(m, s) = s = 
m/3. □ 
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X 



w 



Figure 6: Example of 6-cycle on two weight-2 circulants. 



Lemma 5.9. There is a 6-cycle in case\2\if and only if 

s{p 2 ) = ±2s(p l ) 

Proof. It can be assumed that there is a 6-cycle if and only if if both (cycle) column 
y and column z lie in C 1 and column v lies in C 2 (Figure©. 

Applying Proposition 12.91 4 to cycle column y and cycle column z yields 



(14) 



x-t = ±s] 



(15) 



x — w = =Fs • . 



Applying Proposition 12.91 -2 to cycle column v yields 



(16) 



w-t = ±si 



Since x - t = (x - w) + (w - t), from (O, £[5]> and (fTBT ): 

± S 1 = zps 1 + ± S 2 , 

from which the desired result is obtained: 



S 2 EE 2S 1 



Lemma 5.10. There is a 6-cycle in case\3\if and only if 



S 1 ± S 2 ± S 3 EEE . 



□ 
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Figure 7: Example of 6-cycle on three weight-2 circulants. 

Proof. It can be assumed that there is a 6-cycle if and only if column y lie in C 1 , 
column z lie in C 2 and column v lies in C 3 (Figure IT]). 
Applying Proposition 12.91 2 to cycle column y yields 

(17) x-t = ±s l , 
Applying Proposition 12.91 -2 to cycle column z yields 

(18) t-w = ±s 2 , 
Applying Proposition 12.91 -2 to cycle column v yields 

(19) w-t = ±s 3 , 

Since (x - t) + (w - x) + (t - w) = 0, from £T7]>, (HQ and ( fT9l : 

±s 1 ±s 2 ts 3 = 0, 

from which the desired result is obtained. Note how the sign of one of the separa- 
tions, s 1 in the main theorem, can be fixed. In fact if — s 1 + s 2 — s 3 = then also 
s 1 — s 2 + s 3 = and this is part of the conditions. 

□ 

Lemma 5.11. There is a 6-cycle in case^if and only if 

e(p 2 ) - e(p 2 ) = ±s(p 1 ) ± s(p 3 ) 

Proof. It can be assumed that there is a 6-cycle if and only if simultaneously cycle 
column y lies in C 1 , cycle points (x, z) and (t, v) lie in A 2 and cycle row w lies in 
C 3 (Figure [8]). 

Since cycle column y lies in C 1 , applying Prop. 12.91 2: 

(20) x - t = ±s 1 . 
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Figure 8: Example of 6-cycle on three weight-2 circulants. 



Since cycle row w lies in C 3 , applying Proposition 12.91 3: 

(21) v - z = ±s 3 . 

Since cycle points (x, z) and (t, v) lie in A 2 , applying Lemma l3~4l 

(22) z = x + e 2 , v = t + e 2 . 
Substituting (l22l and (l20l) into (1211) the desired result is obtained. 

□ 

Lemma 5.12. There is a 6-cycle in case\5\if and only if 

Proof. It may be assumed that there is a 6-cycle if and only if, simultaneously, 
cycle column y lies in C 1 , cycle point (x, z) lies in C 1 , cycle point (w, z) lies in 
A 3 , cycle point (w, v) lies in A 4 and cycle point (t, v) lies in A 2 (Figure©. Since 
cycle column y lies in C 1 , applying Proposition 12.91 -2: 

(23) x - t = ±s 1 . 
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Figure 9: Example of 6-cycle on one weight-2 and three weight-3 circulants. 

Since cycle point (x, z) lies in C , applying Proposition 12.91 1: 

(24) z = x + e(p l ). 

Since cycle points (w, v) lies in A 4 , applying Proposition 12.91 1: 

(25) v = w + e(p 4 ) . 

Since cycle points (w, z) lies in A 3 , applying Lemma [3~4l 

(26) z = w + e(p 3 ) . 

Since cycle points (t, v) lies in A 2 , applying Lemma [34] 

(27) v = t + e(p 2 ) . 

The desired result is obtained from (|23i (1241 ). d25l) . d26l and (JSJ. □ 
Lemma 5.13. There is a 6-cycle in case^if and only if 

e( P 2 )-e(p 3 )-e(p A ) + e(p B ) = ±s{p 1 ) 
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Proof. It may be assumed that there is a 6-cycle if and only if, simultaneously, 
cycle column y lies in C 1 , cycle point (x, z) lies in A 2 , cycle point (w, z) lies in 
A , cycle point (w, v) lies in A 5 and cycle point (t, v) lies in A 3 (Figure IT0l>. Since 
cycle column y lies in C 1 , applying Proposition 12.91 -2: 

(28) x - t = ±s 1 . 

Since cycle points (x, z) lies in A 2 , applying Proposition 12.91 1: 

(29) z = x + e(p 2 ) . 

Since cycle points (w, v) lies in A 5 , applying Proposition 12.91 1: 

(30) v = w + e(p 5 ) . 

Since cycle points (w, z) lies in A 4 , applying Lemma [341 

(31) z = w + e(p 4 ) . 

Since cycle points (t, v) lies in A 3 , applying Lemma l3~4l 

(32) v = t + e(p 3 ) . 

The desired result is obtained from ([28]), (|29]>, ([30]), (EB and 021). □ 
Lemma 5.14. There is a 6-cycle in case\7\if and only if 

e(p l) _ e(p 2 ) _ e (p3) + £(p 4) + e(p 5 } _ e(p6) = 

Proof. It can be assumed that there is a 6-cycle if and only if, simultaneously, cycle 
point (x, y) lies in A 1 , cycle point (t, y) lies in A 3 , cycle point (x, z) lies in A 2 , 
cycle point (u>, z) lies in A 5 , cycle point (w, v) lies in A 6 and cycle point (t, v) lies 
in A 4 (Figure ITTb. Since cycle points (x, y) lies in A 1 , applying Proposition 12.91 1: 

(33) y = x + e^ 1 ) . 

Since cycle points (t, y) lies in A 3 , applying Proposition 12.91 1: 

(34) y = t + e(p 3 ) . 

Since cycle points (x, z) lies in A 2 , applying Proposition 12.91 1: 

(35) z = x + e(p 2 ) . 

Since cycle points (w, v) lies in A 6 , applying Proposition 12.91 1: 

(36) v = w + e(p 6 ) . 

Since cycle points (w, z) lies in A 5 , applying Lemma l3~4l 

(37) z = w + e(> 5 ) . 

Since cycle points (t, v) lies in A 4 , applying Lemma [341 

(38) v = t + e(p 4 ) . 

The desired result is obtained from 03]), (O, CH), ([36]), 03) and OH)- □ 
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Figure 10: Example of 6-cycle on one weight-2 and four weight-1 circulants. 
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Figure 11: Example of 6-cycle on six weight-1 circulants. 



44 



It has hence been proved that the conditions listed on the statement considers 
all the possible 6-cycles that can exist on the studied quasi-cyclic matrices. 

□ 

Following two examples. In the first the polynomial p(x) = 1 + x 2 with 
m = 6 for such polynomial s(p) = 2 hence s(p) = m/3 and for condition 1 a 
6-cycle exist. 
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In the second example the first polynomial is p 1 (x) = 1 + x 2 with m = 7 for such 
polynomial s(p l ) = 2 and the second is p 2 (x) = x + x 5 for which s(p 2 ) = 4 
hence s 2 (p) = 2s 1 (p) and for condition 2 a 6-cycle exist. 
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5.3 Conditions for the 8-cycles 

The following lemma reduces the number of cycle configurations to be considered. 

Lemma 5.15. For any matrix M € C ma g^ s.t. at least two weight-2 circulants 
lie in the same row or column then the girth is less or equal to 8. 

Proof. The case of two weight-2 circulants laying in the same d.r is considered, 
the other case is the transpose of such case, when two weight-2 circulants lie in the 
same d.r. The following cycle configuration always arises 

| c 1 - 4 C 2 - 4 | . 

It can be assumed that there is a 8-cycle if columns y and u lie in the first C and 
columns v and z lie in the second C. Cycle column y and cycle column u satisfy 



(39) 



t = ±s l 



(40) 



w 



I = ±s 1 
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Figure 12: Example of 8-cycle on two weight-2 circulants. 



Cycle column v and cycle column z satisfy 

(41) t-w = ±s 2 , 

(42) x - I = ±s 2 . 

Combining equations (l39l >. (l40l>.(l4Tb.(l42T) an equation of the form is 1 ± s 1 ± s 2 ± 
s 2 = is obtained. In particular consider the case s 1 — s 1 + s 2 + s 2 it is equal to 
independently to the values of the two separations. □ 

The previous lemma shows how any H matrix of a LDPC quasi-cyclic code 
may have girth bigger than 8 if and only if, or it contain only weight- 1 circulants or 
the weight-2 circulants are not in the same row or column. Hence the lemma gives 
a practical guideline to the construction of quasi-cyclic codes with high girth. 

The following theorem lists all the possible configurations that may contain 
cycles of length 8. Beside each configuration the conditions on separations and 
exponents under which such cycles exist are given. 
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Theorem 5.16. Let be M G C m , ai p~. The configurations in M that may contain 
a cycles of length exactly 8, are the following^ 

1. 

|C-8|, s{p) = m/A 



C l -5 J 2 -I 

J3_! A 4_] 

e(p 1 )-e(p 2 )-e(p 3 ) + e(p 4 ) = ±2 S (p 1 ), 



C 1 - 4 J 2 - 2 
C 3 - 2 

±s(p 3 ) = 2s(p 1 ) 



C x -3 J 2 -l 
J 3 -l C 4 -3 

e(p l) _ e (p2) _ e(p 3) + e(p 4) = ±s(p l) ± s(p3)j 



e(p 1 ) + e(p 1 ) - e(p 2 



A 1 -2 A 2 -2 
A 3 -2 A 4 -2 

<P 2 ) ~ <P 3 



e(p 3 ) + e(p 4 ) + e(p 4 ) = 0, 



C 1 - 4 J 2 - 1 J 3 - 1 
O A 4 -l A 5 -l 

e(p 2 )-e(p 3 )- e (p 4 ) + e(p 5 ) = ±2 S (p 1 ), 



7. 



C 1 - 3 O J 2 -I 
J 3 - 1 C 4 - 2 J 5 - 1 

e (pl) _ e (p2) _ e(p 3) + e(p 5) = ±s(p l) ± a(p 4) 9 



7 Configurations with two or more weight-2 circulants in the same row or column are not listed 
since they always contain a cycle of at most 8 (Lemma [5.15l >. hence they do not add any new infor- 
mation. 
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A 1 - 2 A 2 - 1 A 3 - 1 
A 4 -2 A 5 -l A 6 -l ' 

e(p l) + e(p l) _ e(p 2) _ e(p 3) _ e( /) _ e( /) + e(p 5) + e(p 6) = , or 
e(p 2) _ e(p 3) + e(p 4) _ e(p 4) _ e(p 5 } + £(p6) = ±s(p l) j or 
e(p 2) _ e( p3) + e(p l) _ e(p l) _ e(p 5) + e( /) = ±s(p 4) . 

C 1 -2 J 2 -2 O 
O J 3 -2 C 4 -2 ' 

C 1 - 2 O J 3 - 1 J 4 - 1 
O C 2 -2 J 5 -l J 6 -l ' 

e(p 3) _ e( p4) _ e(p 5) + e(p 6) = ±s(p l) ± s( p2 }) 

A 4 -l A 2 -l A 3 -l A 4 -l 
A 5 -l A 6 -l A 7 -l A 8 -l ' 

e(p x ) + e(p 2 ) - e(p 3 ) - e(p 4 ) - e(p 5 ) - e(p 6 ) + e(p 7 ) + e(p 8 ) = , or 
tip 1 ) - e(p 2 ) + e(p 3 ) - e(p 4 ) - e(p 5 ) + e(p 6 ) - e{p 7 ) + e(p 8 ) = , or 
e(p l) _ e(p 2 ) _ e(p 3 ) + e (p4) _ e( p5) + e( /) + _ e(p 8) = o . 



C 1 - 3 J 2 - 1 O 
J 3 -l O A 4 -l , 
O A 5 -l A 6 -l 

e(p l) _ e(p 2 ) _ e(p 3 ) + e (p4) + e( p5) _ £(p 6) = ±s(p l )> 



ji _ 2 A 2 - 1 A 3 - 1 
C 4 -2 O O 
O A 5 -l A 6 -l 

e(p 2 )-e(p 3 )-e(p 5 ) + e(p 6 ) = ±s(p 4 ) 
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14. 



A 1 - 2 A 2 - 1 A 3 - 1 
A 4 -l A 5 -l O 
A 6 - 1 O A 7 - 1 

e(p i) + e(p i) _ e(p 2) _ e(p 3) _ e(p 4 } + e(p5) _ e(p6) + e(p 7 } = ; or 

e(p 2) _ e( p3) + e(p 4) _ e(p 5) _ e( /) + e( /) = ±s(p l )> 



75. 



C 1 - 2 O O 
J 2 -I A 3 - 1 O 
J 4 - 1 J 5 - 1 C 6 - 2 

e(p 2) _ e(p 3) _ e(p 4) + e(p 5) = ±s(p l) ± s(p 6 )) 



75. 



A 1 -2 A 2 -l A 3 -l A 4 -l 
A 5 - 1 A 6 - 1 O 
O O A 7 -l A 8 -l 

e(p 2) _ e( p3) _ e(p 4) + e(p 5) + e( /) _ e( /) = ±s(p l )> 



77. 



C 1 -2 J 2 -l J 3 -l O 
O A 4 -l O A 5 -l 
O O A 6 -l A 7 -l 

e(p 2) _ e( p3) _ e(p 4) + e(p 5) + e( /) _ e( /) = ±s(p l )> 



78. 



A 1 -! A 2 -l O O 

A 3 -l O A 4 -l O 

O A 5 -l O A 6 -l 

O O A 7 -l A 8 -l 

e(p l) _ e(p 2) _ e( p3) + e(p 4) + e(p 5) _ e(p 6) _ e( /) + e(p 8) = 0, 



Proof. Proving this theorem is particularly long, the arguments used are similar to 
the ones used to prove the two previous theorems. The detailed proof can be found 
in Appendix. 

□ 

The theorems above listed all the conditions (regarding the exponents and sepa- 
rations) that must be satisfied for a cycle of length less than ten to exist on a generic 



49 



quasi-cyclic H matrix. Checking all the conditions listed in a H matrix can be a 
daunting task. The complexity raises from the fact that for every configuration all 
the possible row and column permutations must be considered. The following sec- 
tions present several classes of quasi-cyclic codes that have some regularity and 
reduce the total number of conditions needed. 

6 A first subclass of codes 

In this section a first family of quasi-cyclic matrices is presented. The structure of 
this class has been developed to achieve high sparsity and to ease implementation 
aspects. This class of codes contain the matrices M in which each d.c. can have 
at most one H and/or one /, where I is the identity matrix. Notation D m>OCj p tl 
denotes such a class. Theorems 13.141 13.151 and 13.161 can be applied to D m ^ a ^ a . 
The use of I sub-matrices reduces the number of variables present in the conditions 
to be check. 

Theorem 1531 15 .7 1 15 . 1 61 are applied to D m ^p a . 

Proposition 6.1. Let M £ D mta ^^. The d.m.'s in M that contain cycles of length 
< 8 and the associate conditions are the following: 

1. Containing A-cycles: 



1.1 



H -4 



s(p) = m/2 



1.2 



H 1 -2 




^P 1 ) = S (P 2 ) 



1.3 



H 1 -! 
I -I 



H 2 -l 



I -I 



eip 1 ) + e(p 2 ) = 



1.4 



H 1 -! 



H 2 -l 



tip 1 ) ~ e(p 2 ) = 



I -I 



I -I 



2. Containing Q-cycles: 



2.1 



H-6 



s(p) = m/3 
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2.2 



2.3 



2.4 



2.5 



2.6 



2.7 



2.8 



2.9 



2.10 



2.11 



|# 1 _4 77 2 - 2|, s(p 2 ) = ±2s(p 1 ) 

Iff 1 _2 77 2 -2 77 3 - 2|, sQo 1 ) ±s(p 2 ) ±s(p 3 ) = 



77 1 - 2 7-2 
O H 2 - 2 



77 1 - 3 H 2 — I 
7-1 7-1 



77 1 - 3 I - 1 
I - 1 77 2 - 1 



#i _ 2 H 2 — 1 H 3 - 1 
O I- 1 7-1 



77 1 - 2 H 2 - 1 7-1 
O I - 1 77 3 - 1 



77 1 - 2 7-1 7-1 
O 77 2 - 1 77 3 - 1 



H 1 -1 H 2 -1 O 
7-1 O H 3 - 1 
O 7-1 7-1 



77 1 - 1 7-1 O 
7-1 O 77 3 - 1 
O H 2 - 1 7-1 



Kp 1 ) = Kp 2 ) 



e (p 1 )- e (p 2 ) = ± s (p 1 ) 



e(p 1 ) + e(p 2 )^±s(p 1 ) 



e(p 2 ) ~ e{p 3 ) = ±s(p 1 ) 



e(p 2 ) + e(p 3 ) = is^ 1 ) 



e(p 2 ) - e(p 3 ) = ±s(p l ) 



e(p l)_ e(p 2) + e(p 3) = 



e (pl) + e (p2) + e (p3) = 



Proof. The theorem can be demonstrated using a case by case analysis of the con- 
figurations presented in previous sessions. □ 
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7 Bresnan Codes 



A variation of the previous class of codes that first appeared in [36 ] and later stud- 
ied by mean of Grobner bases in OBI [371, is investigated. This class is being 
used to reduce hardware complexity of the decoder in 11381 . Such class of codes 
is here referred to with the name "Bresnan codes", they consist of two [a x a] 
square blocks of circulants. Each square block has weight-2 circulants in the main 
diagonal and identities in another diagonal. The structure of the parity check ma- 
trix allows the codes to be (3, 6) -regular LDPC codes and reduces the number of 
conditions that need to be satisfied to have girth 8. 

Definition 7.1 (Bresnan codes). Let a,m be positive integers such that a,m > 4. 
Then 7i m . a denotes the class of the (ma x 2ma) matrices of the form 

H\ ... / 
/ Hi ... 
/ '■■ '■■ : 



H 













h\ 



Hi 




I 

Hi 




I 









I 

















I 

H^ 



where every Hf, with c £ {1, 2} and i £ {1, 2, .., a}, is anmxm binary weight-2 
circulant matrix and I is the Tfi x Tfi identity matrix. Given any matrix H £ ^~im^Q> 
pi denotes the polynomial of H? 

The following fact is obvious. 

Fact 7.2. Let a, m be positive integers such that a,m > 4, then 

H. m ,a — C mjQ , 5 2,l 

The following theorem, presented in [39 ] and later generalized in 041 . gives a 
condition to ensure that H has full rank. 

Proposition 7.3. Let a,m be positive integers such that a, m > 4. For any matrix 
H in H m ,a, let C be the [N, K, d] quasi-cyclic code with parity-check matrix H. 
Suppose 

g cd(i+ n P ix m +i)=i 



Kh<a 



Then N = 2K. 



To further simplify the notation, the sentence "d.s. A is under d.s. B" means 
that both A and B belong to the same d.c. but that the d.r. of A is d.r. ia and the 
d.r. of B is d.r. is, with i& > %b is used. Similarly for "A is at the right of B" and 
other intuitive positional expressions. 

The Bresnan Codes are a subclass of the class presented in section [6] and it is 
possible to apply theorem 16.11 to it. It is clear that any configuration presented in 
the theorems is possible except for 2.8, 2.10-2.14. Moreover for configurations 2.9 
only these special cases can arise: 



52 



1. 







7?-l 



7} +1 -l H) +1 -2 



, j G {1, .., a - 1}, e(pj) + e(p? +1 ) = ±s(p)+i) 



2. 



77) - 1 7T 2 - 2 7? - 1 



j G {1, ..,a - 1}, e(pj) + e(pj +1 ) = ±s(p 2 ) 



H{ - 2 7 1 — 1 77^-1 
77i -1 72-1 



e(pi)+e(p?) = ± S (pi) 



7j - 1 77 2 - 1 



e(p 1 a )+e(p 2 1 )^±s(pl). 



The following important theorem provides a complete characterization for cy- 
cles of length < 8 in H m ,a- A definition is needed in order to present compact 
statements. 

Definition 7.4. Let a > 4 then J a denotes the sub-set of~N 2 such that G J a 
if and only ifi,j G {1, . . . ,a},forj = i + 1 mod (a). 

Proposition 7.5. Let a,m > 4 and M G 7Y miQ . 

• There is a cycle of length 4 in M if and only if at least one of the next 
conditions holds: 

1. m is even and there is 1 < i < a and c G {1, 2} such that 




2. there is 1 < i < a s.t. 

s{p\) = s(p 2 ). 

3. there is G J a s. t. 

e{p})+e{p}) = Q. 

• There is a cycle of length 6 in M if and only if at least one of the next 
conditions holds: 

1. m is divisible by 3 and there is 1 < i < a and c G {1, 2} such that 




53 



2. there is 1 < i < a s.t. 

s{p 2 ) = ±2s{p\), or s{p\) = ±2s{p 2 ), 

3. there is G J a s. t. 

s{Pi) = s{p}) , 

4. there is G J a s. t. 

s{pf) = s(p 2 j) , 

5. there is G J a s. t. 

s{p\) = s(p 2 j) , 

6. there is G J a s. t. 

e(p}) + e(p 2 )^±s(p}), 

7. there is G J a s. t. 

e(p})+e(p 2 ) = ±s(p 2 ), 

8. there is G J a s. t. 

e(pl)+e(p 2 j ) = ±s(p 1 j ), 

9. there is G J a s. t. 

e{p}) + e{p 2 ) = ±s{p 2 ). 

At this stage, it is finally possible to write down the theorem that allows the 
construct of Bresnan Codes with girth at least 8. 

Theorem 7.6. Let a,m > 4 and M G H m ,a- Let g be the girth of the Tanner 
graph of M. Then g > 8 if and only if all the following conditions hold: 

1. for any 1 < i < a and c G {1, 2} 

s( P t) + j, *(J® + f , 

2. for any 1 < i < a, 

s{p1) + s(p]), s(p 2 ) ± ±2s(p}), s(p}) # ±2s(p 2 ), 

3. for any G J a , 

s(pl) sip]), s(p 2 ) s(p 2 ), s(p}) s(p 2 ), 

e(Pi ) + <P)) # 0, 

e(Pi) + <P 2 j) # ±a(pj), <P\) + <P]) # ±«(Pi), 
e(pj) + <P 2 j) ± MP}), e(pj) + e(p 2 ) ± ±s(p 2 ). 
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7.1 Existence of Bresnan Codes 

An exhaustive search has been carried out to show that collections of circulants 
that satisfy Theorem 17 . 61 exist. First some value of m can be excluded because they 
cannot have solutions. 

Proposition 7.7. Let a > 4 and M € 'H m a- There does not exist any M £ TC m ,a 
with m < 10 such that g > 8. 

Proof. The aim of the following proof is to show that with m < 10 the number of 
conditions is too big to allow any solution. The conditions related to the sum of 
exponent in 17.61 will be considered. For sake of simplicity the various separations 
are called s(pj) = a, s{pf) = b, s(pj) = c, s(pj) = d, note also that the conditions 
of Theorem 17 . 61 imply that {a, b, c} are different and that {a, d, c] are different but 
they allow b = d. The second exponent is written as the sum between the first 
and the separation {i.e. e\{2) = e\(l) + a), with this notation it is possible to 
write the condition e\{2) + e|(l) ^ b, as e\(l) + e](l) ^ b - a, The list of all 
conditions related to the sum of exponents in theorem 17^61 using the new notation, 
is the follow: 

ej(l)+e?(l) ^0,a,6,c,d, 

eKl)+ef(l)^-a,-e,-6,-d, 

ej(l) + e|(l) 7^ b — a,c — a,d — a, —2a, —b — a,—c — a, —d — a 
ej(l) + e|(l) a — c, b — c, d — c, —a — c,—b — c, —2c, — d — c 
e\(l) + e| (1) 7^6 — a — c, c — a — c, d — a — c, 

(1) + e|(l) 7^ —2a — c, — b — a — c, —2c — a, — d — a — c 

The conditions are symmetrical for a and c so it can be supposed, without loss 
of generality, that a < c. 

The case with b 7^ d is considered first. Under the condition that {a, 6, c, d} are 
different and considering that the separations are less than m/2, the first two rows 
of the previous list form nine independent conditions that must be satisfied. They 
can be satisfied if and only if m > 9. To prove that there are no solution for m = 10 
it is necessary to find another condition that is independent from the previous nine. 
Condition e}(\) + e){\) ^ b — a is chosen, this new condition is independent from 
the previous nine if b — a is different from {0, a, b, c, d, —a, —b, — c, — d}. The 
following investigate what happen if b — a is equal to any of these. 



b- 


a 


= 


=>b = 


a Not possible 












b- 


a 


= a 


=>b = 


2a Not possible 


b- 


a = 


—a 


=>b = 


Not possible 


6- 


a 


= b 


=> a = 


Not possible 


6- 


a = 


-b 


= > a = 


2b Not possible 


6- 


a 


= c 


=>b = 


a + c 


6- 


a = 


— c 


=> a = 


c + b => a < c Not possible 


b- 


a 


= d 


=>b = 


d + a 


b- 


a = 


-d 


=> a = 


6 + d 



Some cases have been marked impossible because, if they occur, they will break 
conditions related to the separation, in Theorem 17. 61 Such cases cannot give codes 
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with girth higher than 8 so there is no need to study if the new condition is inde- 
pendent. Next it is analyzed what happen in the cases that cannot be discharged, 
considering that with m < 10 {a, b, c, d} < 4. 

b=a+c The conditions a < c and a / c imply that b > 3. 

• If b = 3 then a = 1 and c = 2 but in this case ej(l)+e^(l) 7^ —2c— a is 
an independent condition, in fact the set {0, a, b, c, d, —a, —b, — c, — d} 
became {0, 1, 3, 2, 4, 9, 7, 8, 6} and —2c — a = 5. Note that d must be 
4 to be different from {a, o, c} and d < m/2. 

• If b = 4 then a = 1 and c = 3 (since a ^ 2) but in this case ej(l) + 
e^(l) ^ —2a — c is an independent condition, 

b=d+a The condition a ^ d implies that b > 3. 

• If b = 3 then or a = 1 and d = 2 or a = 1 and d = 2 

- if a = 1 and d = 2 then does not exist c such that c > a, c ^ d, 
c^b and c 7^ 2a, 

- if a = 2 and d = 1 then c = 4 (since c 7^ b) and ej(l) + e^(l) / 
m — d — c is an independent condition. In fact with this values the 
conditions became ej(l) + e|(l)& / {0, 1, 2, 3, 4, 9, 8, 7, 6} but 
m - d - c = 10 -1-4 = 5 

• If 6 = 4 then or a = 1 and d = 3 or a = 3 and d = 1 (since a 7^ d) 

- if a = 1 and d = 3 then c = 2 (since c 7^ b) in this case ej(l) + 
e^(l) 7^ m — d — c is still an independent condition, 

- if a = 3 and d = 1 then does not exist c such that c > a and c 7^ b. 

a=b+d Considering that a < c < m/2 then or a = 2 or a = 3. 

• if a = 2 then 6 = d = 1 that is not possible since b 7^ 2a. 

• if a = 3 then d 7^ 2 since c 7^ 2d so d = 1 and b = 2, in this case 
ej(l) + e|(l) 7^ m — d — c is still an independent condition, 

It has been shown how if d 7^ b it always exist a new independent condition. 
Hence a group of ten conditions that must be satisfied exist, this implied that for 
M G TC m ,a and b 7^ d, it is possible to have 5 > 8 if and only if m > 10. 

Next step is to find a group of ten independent condition in the case b = d. If 
b = d the first two rows of |43]reduce to seven independent conditions so that they 
cannot be satisfied for m < 8. To prove the proposition it is necessary to find three 
more independent conditions. 

The condition e\(l) + e^(l) 7^ b — a is reconsidered when b = d. Following 
the previous reasoning this is an independent condition if b — a is different from 
{0, a, b, c, —a, —b, — c}. 
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The following studies what happen if 6 — a is equal to any of these. 
6 — a = => b = a Not possible 

b — a = a => b = 2a Not possible b — a = —a => 6 = Not possible 

b — a = b => a = Not possible b — a = —b => a = 2b Not possible 

b — a = c => b = a + c b — a = c => a = c + 6 => a < c Not possible 

For the case b = a + c, equivalent to the cases (a = 1,6 = 3,c = 2) and 
(a = 1,6 = 4,c = 3), it is easy to verify that condition —2a — c is an independent 
condition for m < 10. Note how for m = 8 the case a = 1,6 = 4,c = 3 is not 
allowed since 6 = to/ 2. 

Hence e*(l) + e|(l) / 6 — a (or e*(l) + e^(l) ^ to — 2a — c in the case 
6 = a + c) is a new independent condition for to < 10. A set of eight independent 
conditions for m < 10 has been obtained, another two are still needed. The two 
conditions needed must be independent for to = {9, 10} but not necessary for 
to = 8, since it has been proved that for m < 8 there cannot be Bresnan codes 
with girth more than 8. 

Condition e|(l) + e^(l) ^ b — c is considered next, this new condition is inde- 
pendent from the previous eight if d — c is different from {0, a, b, c, —a, —b, —c, 6 — 
a}. Note the case b = a+c for which the eight conditions are {0, a, b, c, —a, —b, — c, —2a— 
c} ( and not {0, a, b, c, —a, —b, —c, b — a}) will be considered, together with other 
particular cases, at the end. 

The following considers what happen if b — c is equal to any of these, consid- 
ering that d = b 

b — c = => d = c Not possible 

b — c = a => b = a + c b — c = —a => c = a + b 

b — c = b => c = Not possible b — c = —b => c = 2d Not possible 
b — c = c => d = 2c Not possible b — c = —c => 6 = Not possible 
6 — c = 6 — a => a = c Not possible 

As mention above, the case 6 = a+c will be treated apart, the only case to consider 
is when c = a + 6 . 

• c = a + 6; to, < 10 implies c < 4 

- if c = 2 then a = 6 that is not allow, 

- if c = 3 then a = 26 or 6 = 2a that is not allow, 

- if c = 4 then the only two possible situations are a = 1,6 = 3 or 
a = 3, 6 = 1. These are not possible if to, = 9 because one separation 
is equal to to/3. In the case m = 10 e\(l) + e^(l) ^ m — a — b — c is 
a new independent condition. 

This proves how e|(l) + ej (1) / d — c is an independent condition for to < 10 
apart for the cases (a = 1,6 = 3,c = 4) and (a = 3,6 = l,c = 4) in which case 
ej(l) + ej(l) + -a-b-cis. 
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At this point it is possible to build a set of nine independent conditions for 
m < 10 . Hence it is proved that for m < 9 there cannot be a code with girth> 8. 
It is now necessary to prove the existence of at least another independent condition 
for m = 10. Condition ej(l) + ej(l) + 6 — a — c is now considered. The 
situations when this new condition is dependent to the condition ej(l) + ej(l) / 
0, a, b, c, —a, —b, — c, b — a, b — c are studied next. In other word situations when 
b — a — c is equal to any of 0, a, b, c, —a, —b, — c, b — a,b — c are considered . The 
special cases (a = 1,6 = 3,c = 4) and (a = 3,6 = l,c = 4) will be considered at 
the end. 

b — a — c = => b = a + c 

b — a — c = a => b = 2a + c b — a — c = —a => b = c Not possible 

b — a — c = b => a = —c Not possible b — a — c = —b => 2b = a + c 

b — a — c = c => b = a + 2c 6 — a — c = — c => 6 = a < c Not possible 

The case b = a + c will be treated a part and it is of no concern at the moment. 
What happen in the cases that cannot be discharged is studied; remembering that 
with m < 10. 

b=2a+c since a < c and b < 5 b is 4. This implies a = 1 and c = 2. This is not 
possible since c = 2b = 2d 

2b=a+c since a < c and b < 5 it is : 

- if b = 2 then a = 1 and c = 3. This is not possible since b = 2a 

- if b = 3 then a = 2 and c = 4. In this case —b — a is a new independent 
condition 

- if b = 4 this is not possible since does not exist a, c s.t. a + c = 8 and 

a < c < 4. 

b=2c+a this is impossible since a < c implies b = 2c + a > 4 

This proves how e\{\) + e|(l) 7^ b — a — c is an independent condition for m < 10 
apart for the case (a = 2,6 = 3,c = 4) in which case ej(l) + e|(l) 7^ — 6 — a is. 

It has so been proved that the set of conditions 
e\{\) + e|(l) 7^ 0, a, 6, c, —a, —6, — c, b — a,b — c,b — a — c 
(or in the case (a = 2,6 = 3,c = 4) the set e l 1 (l)+e|(l) 7^ 0, a, 6, c, —a, —6, — c, 6— 
a, 6 — c, —6 — a) consist on a set of independent inequalities that must be satisfied, 
hence this required that m > 10. 

To prove the proposition it is necessary to consider the four special cases that 
has been neglected : (a = 1,6 = 3,c = 2), (a = 1,6 = 4,c = 3), (a = 
1,6 = 3,c = 4) and (a = 3,6 = l,c = 4) and demonstrate that for any of these 
cases exist three independent conditions to add to the basic set e\{\) + e|(l) 7^ 
0, a, 6, c, —a, —6, — c, when m < 10. 

These special cases do not need to be considered for m = 8 since it has already 
been proved that there is no Bresnan codes with g > 8 in such case. Moreover, in 
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all four cases one of the separation is equal to three but this is not allow if m = 9, 
hence any of this cases cannot happen if m = 9. 

It is only necessary to consider these cases for m = 10; every case is considered 
singularly. 

(a = 1, b = 3, c = 2) It is straight forward to verify that the following set of con- 
ditions are independent {0, a, b, c, —a, —b, — c, —2a — c, —a — b — c,—b — c} 

(a = 1, b = 4, c = 3) It is straight forward to verify that the following set of con- 
ditions are independent {0, a, b, c, —a, —b, —c, —2a — c, —a — b — c, —2a] 

(a = 1, b = 3, c = 4) It is straight forward to verify that the following set of con- 
ditions are independent {0, a, b, c, —a, —b, —c, b — a, —a — c, —2a} 

(a = 3, b = 1, c = 4) It is straight forward to verify that the following set of con- 
ditions are independent {0, a, b, c, —a, —b, —c, b — a, —a — b — c, —b — c} 

It has so been proved that for every m < 10 it is always possible to find a set of m 
independent values from whom the sum of exponents must differ. This implies that 
for m < 10 it is never possible to meet such conditions hence to find a Bresnan 
code with girth 8. □ 

Table[T]shows how many solutions exist for any value of a and m. The percent- 
age given beside some values indicates that due to memory overflow the complete 
search could not be completed. The given number is an estimate based on that 
percentage. Table [2] shows the time taken to find all the possible codes with given 
a and m. It can be seen that the number of Bresnan codes increases with the 
dimension of the circulants (m) and with the dimension of the block (a). 
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Table 1 : Number Bresnan codes for given a and m 
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Table 2: Timing (in seconds) for find all Bresnan codes for given a and m 
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7.2 Performance 

A study of the performances of the Bresnan codes when compared with random 
like code is presented here. 

All simulations presented in the remainder of this thesis have been obtained 
using a maximum limit of 20 iterations. For every SNRdB point a minimum of 10 
block errors has been found. Error bars are plotted on plots in this thesis. The error 
associated with the Frame Error rate is defined as ( fiUTl ): 

(44) FER± = FER x exp (j ^ZZ*™ ) > 

yy n x N error j 

where N is the number simulation trials and N error is the number of frame errors 
recorded. This is to be considered an important aspect of the performance graphs 
often overlooked. 

The simulations show that the codes of the family here developed behave very 
well, even for medium lengths. First the result for a N = 808 code. In Figure [T3l 
a comparison between one of the new codes and three codes taken from known 
good families is provided. All codes are LDPC codes, with rate 1/2 and length 
N = 808. The new code is a regular code with a structure, while the others are 
randomly constructed using some optimization method. To be more precise, these 
three literature codes are: 
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• A Mac Kay code with N = 808, K = 404, d v = 3, d c = 6 obtained with 
the method presented in HD1 ; 

• A randomly constructed code with N = 808, K = 404, d v = 3, d c = 6 this 
code has been obtained by running the optimized edge-placement algorithm 
for five passes (see for example [36 ]); 

• A Richardson irregular code with N = 808, K = 404, this code has been 
created using the optimized weight distributions proposed by Richardson et 
al. MM; 

The performance of the new code is close to that of the random-edge one, with 
the Random code hitting an error floor around 3dbm le — 05. The MACKAY code 
performs worse than the new codes losing 0.5dB at le — 06. The Richardson code 
has particularly poor performance, it is possible that the code constructed is a bad 
representative of the ensemble of which it is part. It is known that for medium and 
shorter length codes the performance of the single codes varies widely from the 
average performance of the ensemble and from the performance of long codes for 
the same construction. On the other hand the situation highlights the problem of 
finding good random codes: creation and simulation of many codes are necessary 
to find a good one. 

Figure [14] shows the performances comparison between ./V = 1048 i? = 1/2 
codes constructed with MacKay method BDl . Random code [41] and a quasi-cyclic 
code of this class. It can be seen that the performance of the quasi-cyclic code is 
close to the one of the random one, it is slightly worse when the code meets a noise 
floor. The performance of the MacKay code in the ldB — 4dB range is worse than 
the quasi-cyclic code. For low values of SNR the quasi-cyclic code performance 
compares well with the random code, and both are superior to the MacKay one. 

Even if the main interest of this thesis is on medium length codes it is interest- 
ing to evaluate the construction presented for longer codes to have an idea of how 
such family performs in such case, for this reason two code with N 10, 000, are 
presented. In Figure [l5]a Bresnan code is compared with a MacKay code, it can be 
seen how the MacKay code has a much better defined waterfall region that make 
the code outperform the Bresnan code for high SNR, this could be due to low min- 
imum distance or to the presence of trapping sets. The result is in line with what 
is a know fact: analytic codes do not perform as well as random codes for high 
length. Moreover it proves that the MacKay codes for this length start to behave as 
it is theoretically proved for infinite length. 

Finally it is important to consider that for the concentration theorem (9l it is 
more likely that a code taken at random from the ensemble performs as well as 
the average performance of the ensemble for this value of N. It is important to 
observe that the construction of good random codes, such as the three confronted 
here, requires a consistent amount of computation, while the new codes can be 
obtained immediately, thanks to their simple structure. 
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Decoding Performance (FER) 
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Figure 13: Block error rate performance for N = 808 Bresnan code compared with 
three random like codes of the same length obtained from MacKay, Richardson and 
random permutations methods. 
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Figure 14: Block error rate performance for N = 1048 Bresnan code com- 
pared with the a code obtained with the random permutation method and one with 
MacKay methods. 
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Decoding Performance (FER) 
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Figure 15: Block error rate performance for N 
code and a MacKay code are compared. 



10, 000. A Bresnan quasi-cyclic 



8 Extension of the Bresnan Codes 

In this section an approach to construct Bresnan-like codes with rate higher than 
1/2 is presented that maintain girth 8. 

The Bresnan codes are formed by two square blocks of circulants, the codes 
developed in this section are formed by three blocks. A rate of 2/3 is achievable 
in such way. The method is expandable to more blocks and any rate of the type 
(I — I) /I can be achieved. 

Definition 8.1 (Rate 2/3 Bresnan codes). Let a, m be positive integers such that 
a > 4, m > 3. M mjQ , denotes the class of the {ma x 3m«) matrices of the form 



H 



Hi 
I 








H\ 















Hi 




I 



I 

Hi 









I 







I 



Hi 



Hi 





Hi 



I 













where every Hf, with c € {1, 2} and i € {1, 2, .., a}, is anmxm binary weight-2 
circulant matrix and I is the m x m identity matrix. 

To simplify the discussion the three blocks that form the H matrix are called 
B\, f?2 and B% and H = [U1I.B2I.B3]. The following theorem lists all the condi- 
tions that must hold for an H matrix of this type to have girth exactly eight. 
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Theorem 8.2. Let m > 3, a > 4 and M £ H mja ,. Lef g be the girth of the Tanner 
graph of M. Then g = 8 if and only if all the following conditions hold: 

1. for any 1 < i < a and c £ {1, 2, 3} 



2. /or any 1 < i < a and c,d£ {1, 2, 3} and c ^ d, 

s(pt) + s(pf), s(pt) £ ±2s(pf), 

3. for any 1 < i < a, 

s(pj) ± s(p 2 ) ± s(pf) # 0, 

4. for any £ J a , 

4.1 e(pj) + e(p|) ^ 0, 

4.2 e(p*) + e(p?) # {±*(p*),±s(p?)}, 

4.5 e(pj) + e(p 2 ) # {± S (p}),± S ( Pl 2 ),± S (pf),± S (p|)}. 

5. for any 1 < i < a, 

s(Pi) # {«(Pi+l)'' S (^+l)' S (^+l)'' S (^+2)} ; 

Kft 2 ) # {s(Ph-i),s(p? +2 )}, 

«(Pi)#{*(Pi+2).«(^+l)} S 

6. for any 1 < i < a, 

6.1 e (pj) + e(pj +1 ) + e(pf +2 ) ^ 

6.2 e ( p 2 )+e(p 2 +i) _ e(p 3 +i)#0 

6.5 6(pi)-6(p 2 ) + 6(pf +1 )^0 
6.4 e (pi)- e ( p 2 +2 ) + e(p 3 +2) ^ 

Proof. It is necessary to prove that the conditions listed in the statement cover all 
the possible cycle configurations existing for the construction presented. To prove 
this it is necessary to show that the configurations, regarding cycles of length 4 and 
6 in theorem 16.11 or are not possible for this construction, or are associated with 
conditions listed in the statement. 

• Configurations 1 . 1 and 2. 1 in theorem 16.11 are evidently associated with the 
conditions in point 1 . 

• Configurations 1.2 and 2.2 are covered by the conditions in point 2. 

• Point 3 considers configuration 2.3 in theorem |6~T1 
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Configuration 1.3 is considered by condition in point 4.1. Note that due to 
the position of the diagonal of identities in B3 such configuration can exist 
only between B\ and B2. 

Configuration 1.4 cannot exist in the proposed class thanks to the position of 
the identities diagonals. There cannot be any couple of columns that have H 
and I sub-matrices in the same two rows. 

Configuration 2.4 in considered by the conditions in point 5. In fact every 
H sub-matrix has three / sub-matrices lying in the same row and one in the 
same column. Hence for every H sub-matrix six possible cycles of the type 
2.4 can exist, three formed by the three identities in the same row and three 
formed by the identity in the same column and the three H sub-matrices 
lying in the row of such /. 

The conditions for cycles to existing in such configuration are: 
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Eliminating from the list above the conditions that are equivalent the set of 
conditions in point 5 is obtained. 

Configurations 2.5 cannot exist for the same reason discussed for configura- 
tion 1.4, the same is true for configurations 2.7 and 2.9 of theorem |6~T1 

Cycles from configuration 2.6 can exist only between sub-matrices in B\ 
and sub-matrices in B2 . In fact the position of the diagonal of identities in 
B3 does not allow any of the columns in the third block to overlap in two 
positions with any of the other columns. Conditions in point 4.2 evidently 
cover all the possible cycles of this type. 



when the 



columns are in B\ and B2. The remaining column 



For the same reason of the previous point, configuration 2.8 can exist only 

H 2 -2 I 

can be in any of the three blocks, the relative cycles are avoided by the con- 
ditions in point 4.3. 



Cycles from configurations 2.10 and 2.11 are associated to conditions in 
point 6. To prove this it is necessary to show that the conditions listed in point 
6 consider all possible cycles that can arise from the two configurations. 

Both configurations require that the distance from an H sub-matrix and the 
identity in the same column is the sum of the distances from the other two 
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H sub-matrices and relative identities. It is evident from the position of the 
identities diagonals that such case is possible only if one column lies in B% 
and the other two columns do not. If the other columns both lie in B\ it is 
only a matter of connecting the sub-matrices to see that the resulting cycle 
is part of configuration 2.11 and it is not allowed by the condition 6.1. If 
the other columns both lie in B2 it is only a matter of connecting the sub- 
matrices to see that the resulting cycle is part of configuration 2.10 and it 
is not allowed by the condition 6.2. If one of the remaining column lies 
in B\ and the other in B2 there are two possibilities to form a cycle, both 
corresponding to configuration 2. 10. The fist case has the H sub-matrices 
in B\ and B2 lying in the same d.r., the second case has the H sub-matrices 
in B2 and B3 lying in the same d.r. The resulting cycles are avoided by 
conditions 6.3 and 6.4 respectively. 

Note that if the two H sub-matrices in B\ and B% lie in the same d.r. is not 
possible to form a cycle due to the position of the identities diagonals 

It has been proved that all the possible cycles that can exist in the H matrix for this 
class of codes is avoided by one of the condition listed. □ 

The performance of such class of codes is presented next. Figure [16] shows the 
comparison between a code from the class of quasi-cyclic codes presented versus 
a random code and a Mac Kay code obtained from [42], the codes have length 
N ?s 200. The codes have similar length and rate, the MacKay code is slightly 
longer but has higher rate R k, 0.7. It can be seen how the quasi-cyclic code and 
MacKay code perform closely and slightly better than the random code. 

The comparison of the same families of codes for N ks 1 , 000 is presented in 
Figure [T7] The codes have similar length and rate, the MacKay code is slightly 
longer but has higher rate R m 0.85. For such length the proposed quasi-cyclic 
code outperforms both the random code and the MacKay code. In particular the 
random code seems to hit an early error floor and the MacKay code perform few 
tens of dB worse for all the SNR range considered. 

Finally Figure [18] presents the comparison of a quasi-cyclic code with a ran- 
dom code for N = 2457. Unfortunately it has not been possible to obtain a 
MacKay code with such rate and length due to failure of the search carried out. 
The random code performs slightly better than thequasi-cyclic code for low SNR 
but it is affected by early error floor that makes its performance poor at higher 
SNR. In contrast the quasi-cyclic code does not show the presence of any early 
error floor. 

It can be concluded that the construction proposed is a valid alternative to ran- 
dom like codes for high rate. The quasi-cyclic codes, for the dimensions con- 
sidered, perform at least as well as other constructions but offer the advantage of 
reduced complexity and minimize memory requirement. Those are important as- 
pects for many applications where high rate codes are used. 
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Decoding Performance (FER) 
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Figure 16: Block error rate performances comparison of a regular quasi-cyclic 
codes with R = 2/3 versus a MacKay codes and a random permutation codes of 
similar rate and length N m 200 
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Figure 17: Block error rate performances comparison of a regular quasi-cyclic 
codes versus a random constructed code and a MacKay code, N ss 1,000. All 
codes have similar rate. 
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Decoding Performance (FER) 
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Figure 18: Block error rate performances comparison of a regular quasi-cyclic 
codes versus a random constructed code with same rate and length ]V«2, 000 

9 Two families of quasi-cyclic codes with girth 10 

In this section two families of quasi-cyclic LDPC codes with girth at least 10 and 
rate 1/2 are presented. The first construction is a (2, 4)-regular codes, the second 
a (3, 6) -regular codes. The constructions are based on the idea of using blocks of 
circulants to compose the H matrix. The circulants lie in the diagonals of the 
block, this help to detect the possible configurations arising and can be exploited 
to reduce hardware complexity of the decoding algorithm (see EH). 

9.1 Regular (2, 4) 

This family of codes has been studied first because its structure drastically reduce 
the number of conditions necessary to ensure high girth. 

Definition 9.1 ( (2, 4)-regular quasi-cyclic codes girth 10). Let a, m be positive 
integers such that a > 4, m > 3. Notation C m . a denotes the class of the (ma x 
2ma) matrices of the form 

u = [b 1 \b 2 ] 

where 



H 
















#i 












































H a 
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and 



Bo 



I 
I 

J 2 o 







. . . Jo 
... Ji 



Ja-1 








I 



Where every Hi, with i £ {0, 1, .., a — 1}, is an m x m binary weight-2 circulant 
matrix, Jj are binary weight-1 circulant matrices and I is the m x m identity 
matrix. 

Reducing the theorems studied in section \5\ for this class of matrices is a long 
task. This can be eased by considering that some of the cycle configurations listed 
in theorems l5.3H5.7|[5.16l can exist only if the distance from the main diagonal and 
the diagonal of J's in B 2 is a particular value. Such distance is called S. 

The following theorem lists all the conditions that must hold for an H matrix 
of this class to have girth at least 10. 

Theorem 9.2. Let m > 3, a > 4 and H € C 

ra,a- Let 9 ^ e me girth of the Tanner 
graph of H . Then g > 10 if and only if all the following conditions hold for any 
0<i<a-L 



*(Pi)^f, *(Pi)^f, s( Pi )^^, 



s(Pi) + s(p i+ s), 



3. 



or 5 ^ a/2, or < 



or 5 ^ ±a/3, or 
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5. 



or 5 ^ ±a/4, or e(p*) + e(p i+s ) + eOi+25) + e(Pi+3s) + e(j>i+4<5) # 



Where the index operations are modulo alpha. 

Proof. It is necessary to prove that the conditions listed cover all the possible cycle 
configurations existing in such class of codes. To prove this it is necessary to show 
that the configurations listed in theorems 15.3115.71 and 15 . 1 6 1 are not possible for this 
consuuction or are associated with a condition listed in the statement. First the 
configurations that can have 4-cycles, listed in theorem [531 are considered. 

• Configuration Q] is associated with the first condition in point 1. 

• Configuration [2] cannot exist because there cannot be two H sub-matrices in 
the same row or column. 

• Configuration [3] can exist in this construction only if 5 = a/2. In any other 
case there cannot be four sub-matrices lying in two columns and two rows. 
If 5 = a/2 condition 3.1 assures that the cycle cannot exist. 

The configurations that can have 6-cycles, listed in theorem 15.71 are considered 
next. 

• Configuration Q] is associated with the second condition in point 1. 

• Configurations [2] and [3] cannot exist because there cannot be two H sub- 
matrices in the same row or column. 

• Configurations [4] and \5\ cannot exist because for this construction every H 
sub-matrix lies alone in a column. 

• For the same reason discusses for configuration 3 of theorem [531 configura- 
tion[6]can exist in this construction only if S = a/2. If 5 = a/2 conditions 
3.2 and 3.3 assure that the cycle cannot exist. 

• Configurations [7]requires that the distance from two sub-matrices in the same 
column is the sum of the two other such distances for the remaining d.c.'s in 
the configuration. It is evident that all the columns must lie in B2 since in B\ 
there is only one sub-matrix in each column. Moreover such situation exists 
only if 5 = ±a/3. If S = ±a/3 conditions 4.1 assures that the cycle cannot 
exist. To determine the sign of the terms in the condition it is sufficient to 
consider that there cannot be two J's in the same column. Hence the J's 
take alternate positions in the cycle configuration hence they all have the 
same sign. 
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The configurations that can have 8-cycles, listed in theorem 15.161 are considered 
next. 

• Configuration Q] is associated with the third condition in point 1. 

• Configurations I2ll3ll4l 171 [l2"land[T5lcannot exist because for this construction 
every H sub-matrix lies alone in a column. 

• For the same reason discussed for configuration 3 of theorem [531 configura- 
tion [5j can exist in this construction only if 5 = a/2. If <5 = a/2 conditions 
3.4 assures that the cycles cannot exist. 

• For the same reason of previous point, configuration [6] can exist in this con- 
struction only if S = a/2. If 5 = a/2 conditions 3.5 and 3.6 assure that the 
cycles cannot exist. 

• Configuration [8] cannot exist because the structure of the H does not allow 
to have 6 sub-matrices lying in two rows and three columns. 

• Configuration [9] is evidently associated with condition 2. 

• Configuration [TOlcan exist only with 6 = a/2 for the same reason discussed 
in previous cases. If 5 = a/2 conditions 3.7 assures that the cycles cannot 
exist. 

• Configuration QT] cannot exist because the structure of the H does not allow 
to have 8 sub-matrices lying in two rows and four columns. 

• Configuration [13] cannot exist because it implies the existence of three non 
zero sub-matrices in the same column or the existence of a column with an 
H and one identity, both cases are not possible due to the particular structure 
of the H matrix. 

• Configuration [14] cannot exist because it implies to have three non zero sub- 
matrices in the same column that is not possible due to the particular struc- 
ture of the H matrix. 

• Configuration [16] cannot exist because it implies to have four non zero sub- 
matrices in the same column or row that is not possible due to the particular 
structure of the H matrix. 

• Configuration [17] requires that the distance from an two sub-matrices in the 
same column is the sum of the two other such distances for the remaining 
d.c.'s in the configuration. It is evident that the columns not containing H 
must lie in Bi- Moreover such situation exist only if 5 = ±a/3. If 5 = 
±a/3 conditions 4.2 assures that the cycle cannot exist. The same reasoning 
applied to configuration [7] of theorem l5"77"l can be used to explain the signs of 
the terms of the condition. 
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• Configurations [T8]requires that the distance from the two sub-matrices in the 
same column is the sum of the three other such distances for the remaining 
d.c.'s in the configuration. It is evident that all the columns must lie in Bi- 
Moreover such situation exist only if 5 = dba/4. If 5 = ±a/4 conditions 5 
assures that the cycles cannot exist. 

It has been proved that all the possible cycles that can exist in the H matrix for this 
class of codes is avoided by one of the condition listed. □ 

Aside from the degenerative cases where 5 = a/2, 5 = a/3 and 5 = a/4 the 
conditions to check are extremely simple. Even in the degenerative cases the list 
of conditions is quite short. Thanks to this, finding exponents and conditions that 
allows girth 10 is an easy task. Following the performances of three of such codes 
are presented. The codes have been obtained avoiding the degenerative cases of 5 
and imposing the conditions on the separations. 

In Figure [19] the block error ratio versus noise performance graph of a small 
code (N = 810) from this family can be seen. The code is compared with a 
random code. Both codes have same length and are (2, 4)-regular. The random 
codes has been obtained with a optimal permutations method and has a girth of 
14 and an average girth of 15.7, the quasi-cyclic codes has girth 10. In the same 
figure the quasi-cyclic code is compared also with a more advanced Progressive 
Edge Grow (PEG) code from fi3l l44l . The code is (2, 7)-regular. It can be seen 
how this code performs as well as the random code and compares evenly with a 
carefully constructed code, this is a good result for such class of codes. 
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Figure 19: Block error rate performances comparison of a quasi-cyclic code com- 
pared with a random code and a PEG code with N rs 800 

Figure [20]presents a comparison of a quasi-cyclic code from this family versus 
a random code in the case of longer codes. Both the codes have N w 6, 500. It 
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Figure 20: Block error rate performances comparison of a quasi-cyclic code com- 
pared with a random code. Both codes are (2, 4)-regular with N sj 6, 500 

can be seen how the performance of the quasi-cyclic code with girth 10 do not hit 
an error floor as soon as the random code. The random code has average girth of 
20.37 but girth of only 6. 

Regular (2, 4)-LDPC codes are not often used in practice due to their poor per- 
formances. As an example of this a comparison between a (3, 6)-regular Gallagher 
code and the PEG code is presented in Figure [2]]lt can be seen how the (3, 6)- reg- 
ular code has up to 5db gain when compared with (2, 4)-regular codes both random 
and quasi-cyclic . For this reason the interest is on build (3, 6)-regular quasi-cyclic 
LDPC codes with girth ten. 



The construction presented here is the result of few considerations. Lemma [5. 151 
shows how any quasi-cyclic LDPC code with girth at least 10 cannot exist if two 
weight-2 circulants lie in the same row or column. For this reason all weight-2 
circulants are placed in the main diagonal of B\ . To keep the column weight equal 
to three it is necessary to insert another weight- 1 circulant for each column and this 
is done inserting another diagonal of identities. The second block B2 cannot con- 
tain weight two circulants hence any row and column must contain three weight- 1 
circulants to form a (3, 6)-regular code. Aiming to maintain some regularity to 
allow reduced hardware complexity of the decoding algorithm, the circulant matri- 
ces have been disposed in three diagonals. To simplify the number and complexity 
of the checks necessary to guarantee girth 10 identity matrices, as special case of 
weight- 1 circulant matrices, have been used where possible.The resulting family 
of codes is the following. 



9.2 Regular (3, 6) 
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Figure 21: Block error rate performances comparison of a (2, 4)-regular quasi- 
cyclic codes versus a (3, 6) -regular code N ss 800 



Definition 9.3 ((3, 6)-regular quasi-cyclic code). Let a, m be positive integers 
such that a > 4, m > 3. Notation Ti) m ,a denotes the class of the (ma x ma) 
matrices of the form 

H = [Si | B 2 ) 



where 
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and 



J ... I 0...0 I 
OJi 0... I 0... I 



I 



I 



I ... I ... J a _i _ 

where every Hi, with i € {0, 1, .., a— 1}, is anmxm binary weight-2 circulant 
matrix and I is the m x m identity matrix, and Ji are weight-1 circulant matrices. 

Some of the cycle configurations listed in theorems 15.3115.71 and 15 . 1 6 1 exist only 
if the distance from the main diagonal and the diagonals of identities in B2 take 
particular values. Such distances are called 62,63, and for simplicity the notation 
£3 — 62 is called 632- Such values can be seen in Figure |22] Moreover the position 
of the diagonal of identities in B\ is called 6\. For the construction presented here 
6\ = a — 1. 

Note that it is not possible to substitute the diagonal of J's with identities. In 
fact in such case 6-cycle configurations will always arise independently from the 
positions of the diagonals. The existence of this configurations carry a parallel 
with the existence of 6-cycles on a weight-3 circulant matrix ll33l . A graphical 
representation of such case is given in Figure [22] It can be seen that the particular 
6-cycle exists if 62 + 632 = 63 but 62 + 63 — 62 = 63 hence the condition is always 
satisfied and the cycle exists independently from the values of the deltas. If all the 
matrix involved in the cycle are identities it is not possible to avoid such cycle. 
Setting the main diagonal to contain J's, and not identities, it is possible to avoid 
such cycle by imposing proper conditions on the exponents. 

A new class of matrices is presented next. The matrices are part of the previous 
class but some limitations on the position of the diagonals of identities are added. 
Such limitations on the values of deltas are not strictly necessary but help to reduce 
the number of possible configurations, and so the complexity of the conditions to 
check. 

Definition 9.4 ((3, 6) -regular quasi-cyclic code girth 10). Let E m)Q denotes the 
class of the (ma x ma) matrices in B miQ that satisfy the following conditions on 
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Figure 22: Example of configuration existing independently from the values of the 
deltas 

the distances of the diagonals. 
(45) 

1 5 X ^ a/2, S x / ±a/3, S x + ±a/4 where x G {2, 3, 32} 

2 5 X ^ where x, y G {1, 2, 3, 32} a«<i x ^ y 

3 ^ ± <5y ± 5 2 7^ 0, where x,y,z G {2, 3, 32} 

4 ± 5y ± 5 Z ± 5^ ^ 0, where x, y,z,w G {1, 2, 3, 32} 

Once those conditions (on the deltas) are verified it is possible to apply the 
theorems studied in section |5]to obtain the set of conditions (on the circulants) that 
assure girth at least 10. 

Theorem 9.5. Let m > 3, a > 4 and H G E mjQ ,. Let g be the girth of the Tanner 
graph of H . Then g > 10 if and only if all the following conditions hold for any 
0<i<a-l: 

1. 

sipi)*% s( Pi )^f, s(p t )^f, 

2. 



11 



2. 



2. 



2. 



3 



2 



1 



s(Pi) # s(pi+5i)> 
s(pi) ^ ±2s{p i+Sl ), 
s(Pi+<5i) # ±2s(pi), 



3. 



s(Pi) # s (Pi+<jJ> w/ie/B x e {2, 3, 32} 



4. 



4.1 e( Pi ) J -e( Pi±Ss2 ) J #0, 

4.2 e(pi) J - e(p i+± 5 32 ) J ^ ±s(pi), 

5. 

5.7 e( Pi )"-e(^ 32 )^0, 

5.2 e( K ) H - e{ Pi -s x ) H + e(p i+1 ) J - e{ Pi ) J £ where x E {2, 3} 



Proof. It is necessary to prove that the conditions listed cover all the possible cycle 
configurations existing in such class of codes. To prove this it is necessary to show 
that the configurations listed in theorems 15.31 15.71 and 15.161 or are not possible for 
this consuuction, or are associated with a condition listed in the statement. First 
the configurations that can have 4-cycles, listed in theorem 1531 are considered. 

• Configuration Q] is associated with the first condition in point 1. 

• Configuration [2] cannot exist because there cannot be two H sub-matrices in 
the same row or column. 

• Configuration [3] cannot exist in this construction and with the conditions 
imposed on the deltas. In fact such configuration can exist only if there 
are four sub-matrices lying in two column and two rows. This can happen 
only if two deltas are equals modulo alpha, or if one of the delta is equal 
to a/2. These cases are not allowed by conditions 1 and 2 in definition 19.41 
hence the cycle configuration cannot exist. 

The configurations that can have 6-cycles, listed in theorem 15.71 are considered 
next. 

• Configuration Q] is associated with the second condition in point 1. 

• Configurations [2] and [3] cannot exist because there cannot be two H sub- 
matrices in the same row or column. 

• Configuration [4] contains two C sub-matrices in the two d.c. hence both the 
columns must be in the B\. Every H sub-matrix can form a cycle only with 
the row above and below itself. All this cases are considered by condition 
2.1. The remaining sub-matrix in the configuration is an identity hence the 
two exponent terms are not present in the final condition. 



78 



• Configurations [5J and [6] cannot exist for the same reason discussed for con- 
figuration [3] of theorem [531 

• Configurations |7]requires that the distance from two sub-matrices in the same 
column is the sum of the two other such distances for the remaining d.c.'s in 
the configuration. This is not allowed by condition 3 in definition 19.41 with 
the exception of the case depicted in Figure [22] Such cycles are avoided by 
condition 4.1. 

The configurations that can have 8-cycles, listed in theorem 15.161 are considered 
next. 

• Configuration [TJ is associated with the third condition in point 1. 

• Configurations 0glEl|gEZlISlIISlIIIlII2lIIlII3andITg cannot exist for the 
same reason discussed for configuration [3] of theorem [5731 . 

• The same reasoning done for configuration 0] of theorem 15771 applies to con- 
figuration [3] and all the possible resulting cycles are avoided by conditions 
2.2 and 2.3. 

• Configuration [9] contains two C sub-matrices in the two d.c. hence both the 
columns must be in B\. The remaining column contains two J sub-matrices 
hence it must lie B?,- The sub-matrices in the column in B2 are spaced at 
distance 5 X with x G {2, 3, 32} hence all the possible cycles are considered 
by condition 3. 

• Configuration [12] requires that the distance from two sub-matrices in the 
same column is the sum of the two other such distances for the remain- 
ing d.c.'s in the configuration. This is not allowed by condition 3 in defini- 
tion [974] Note that the presence of a weight-2 circulant forces one d.c. to be 
in B\ hence the case depicted in Figure [22]does not exist here. 

• Following the same reasoning done for configuration [7]of theorem 1577] Con- 
figuration [T7] exists only if the three columns lying in B2 have the configu- 
ration depicted in Figure [22] and the last column lies in Bl. Such cycles are 
avoided by condition 4.2. 

• Configuration [T8l is associated to conditions in point 5. To prove this it is 
necessary to show that the conditions listed in point 5 consider all possible 
cycles that can arise from the configuration. 

The configuration requires that the distance of two sub-matrices in the same 
column is the sum of the three other such distances for the remaining d.c.'s in 
the configuration. Such configurations are not possible thanks to condition 4 
in definition 19.41 A particular case remains to be considered. If two columns 
of the cycle configuration lie in B\ and two in B2 then the existence of this 
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configurations carry a parallel with the existence of 8-cycles on two weight- 
2 circulants in the same row or column (Figure [121 . It has been proved 
in lemma 15.151 that such cycle exists independently from the value of the 
separations. In this case the delta's take the part of the separation, hence the 
cycle configuration exists independently from the value of delta's. But Bi 
has three diagonals so three different cycles can exist, and such cycles are 
avoided by the three conditions in point 5. 

To better understand this passage consider that the two cycle columns in B\ 
contain an identity and an H sub-matrix each, hence exponents of two H 
sub-matrices must be part of the condition. From Figure [12] it can be seen 
that the cycle exists only if both the cycle columns in B2 contain two sub- 
matrices appertaining to the same two diagonals. There are three possibili- 
ties: one sub-matrix is in the second diagonal and one in the main diagonal, 
one sub-matrix is in the third diagonal and one in the main diagonal and fi- 
nally one sub-matrices in the second diagonal and one in the third diagonal. 
The first two cases lead to condition 5.2 the third to condition 5.1. 

It has been proved that all the possible cycles that can exist in the H matrix for this 
class of codes is avoided by one of the condition listed. □ 

Following the performances of codes of this family for different lengths are 
presented. 

Figure [23] shows the performance of a medium length quasi-cyclic code when 
compared with a random generated code. Both codes are (3, 6)-regular codes of 
dimension N = 2, 294. It can be seen how the quasi-cyclic codes from this family 
performs really close to the random code, but it seems that the quasi-cyclic codes 
hit an error floor at le — 09 BER. 

Even if out of the boundaries of the range of codes that are the main focus 
of this thesis the comparison of quasi-cyclic codes with random codes using 
longer codes is presented next for completeness. Figure [24] and [25] present such 
a comparison. Figure [24] shows the performances of two (3, 6)-regular codes with 
N R3 6, 500 and Figure [25] shows the performances of two (3, 6)-regular codes 
with N « 19, 500 It can be seen the quasi-cyclic codes are quite close to the ran- 
dom codes. With N m 6, 500 only O.ldb, at le — 09, divides the two codes, with 
N 19, 500 the difference increased but remain under the 0.5db, at le — 05, that 
is a considerable achievement for algebraic codes of this length. 

10 Conclusions anf future work 

In this contribution a study on the presence of cycles on the quasi-cyclic LDPC 
codes has been presented. The aim of the study is to identify the configurations 
of circulant matrices that allow cycles of a given length to exist. Once all such 
configurations are known it is possible to associate to each of them a condition, 
on the exponents and separations of the circulants involved in the cycle, that must 
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Figure 23: Block error rate performances of a girth 10 quasi-cyclic code compared 
with a random generated codes. Both codes are (3, 6)-regular with N = 2, 294 
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Figure 24: Block error rate performances comparison of a quasi-cyclic code com- 
pared with a random code. Both codes are (3, 6)-regular with N 6, 500 
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Figure 25: Block error rate performances comparison of a quasi-cyclic code com- 
pared with a random code. Both codes are (3, 6)-regular with N 19, 500 

be satisfied for the cycle to exist. The collection of these conditions gives a way 
to construct quasi-cyclic codes without undesired cycle by ensuring that all the 
conditions are not satisfied. The chapter started with a novel definition of cycles in 
a matrix, opposed to the normal definition that is given on graph. Such a definition 
is used to prove the "isolation" (lemma [3 .101 ). the lemma is extensively used in the 
search of which configurations are valid and can contain cycles. The configurations 
that allows cycle of length less then 10 are found first for the generic case of a 
matrix that can be partitioned in sub-matrices, then in the case of the matrices for 
quasi-cyclic codes. Our approach is not restricted to circulant matrices of weight 
one as is commonly done in literature but considers also circulant matrix of weight- 
2 (weight-3 circulant matrices are of no interest for this search because they have 
internal 6-cycles). The set of conditions for the generic case is found and then 
the conditions are applied for various particular constructions. The first class of 
quasi-cyclic codes studied are the Bresnan Codes 11361 . several aspects of these 
codes are studied and their performances presented. Such class is then extended 
for higher rate codes. Finally two new structures and the relative set of conditions 
that guarantee girth at least 10 are developed and their performances presented. 
The contribution of this work is not only the four classes of quasi-cyclic codes 
with high girth but mainly the creation of the list of conditions that every quasi- 
cyclic codes (of any possible construction method) must satisfy to have a given 
girth. This not only covers all the existing cases but can be used when constructing 
new structures for quasi-cyclic codes to easily identify where the cycles can exist 
and how to avoid them. 

The methodology presented in this thesis could be used to expand the under- 
standing of quasi-cyclic codes. It would be possible to continue the search of 
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codes with high girth and look for the conditions that guarantee girth twelve but, 
in the authors opinion, such work would be long, tedious and would carry little 
benefit. In fact, other parameters of the code, such as the minimum distance, the 
diameter, the minimum stopping set and others, will became dominant and prevent 
improvement in the performance of the code. Of more interest would be to study if 
the methodology presented could be applied to the problem of finding codes with 
good diameter value. It could be possible to develop a set of conditions that en- 
sures the diameter to be a certain value. The study of how these rules intersect with 
the girth conditions could lead to classes of codes that can achieve given girth and 
diameter and offer a good balance between the two. 
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Appendix 

Proof of theorem 13.161 



Theorem 10.1 (13.16b . Let B e M m ,a,/3,r 

The only possible 8-cycle configurations are as follows (l)-presentations): 



1. (1,1), 

2. (1,2), 

3. (1,3), 

4. (1,4), 

5. (2,2), 



6. (2,3), 



7- (2,4), 



|8| 



2.1 6 2 , 2.2 4 4 , 



4 2 2 



2 2 2 2 



6.1 



6.4 



5.1 



5.4 



4 1 1 

1 1 

2 2 2 

2 
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2 



6.2 



6.5 



3 2 1 
1 1 



1 1 
1 1 



6.3 



6.6 



3 10 
112 

2 2 
2 2 



7.1 



2 2 11 
11 



7.2 



2 110 
112 



7.3 



1111 
1111 



(3,3), 
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1 




2 
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9. (3,4), 
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1 


1 







1 


1 


1 


1 


9.1 





1 





1 


, 9.2, 


1 


1 
















1 


1 










1 


1 



10. (4,4), 





1 








1 





1 








1 





1 








1 


1 



Proof. Lemma [331 proves that these are all the possible dimensions that 8-cycle 
configurations can have. It is now necessary to prove that they are also all and only 
valid configurations. 

Following from Theorem 13 .121 the only possible weights vectors for a 8-cycles 

are: 

[8], 

[6,2],J&^lf, 

M)^, [5,1,1,1], 

[4, 4] , [4, 2, 2] , J^<lf, [4, 1, 1, 1, 1] , 

J^f , [3, 3, 1, 1] , JpM>^f , [3, 2, 1, 1, 1] ,[3,1,1,1,1,1], 

[2, 2, 2, 2] , J2^3<lf , [2, 2, 1, 1, 1, 1] , [2, 1, 1, 1, 1, 1, 1] , 

[1,1,1,1,1,1,1,1], 

Configuration[l]A type (1,1) configuration can be generated only by weights vec- 
tor [8] and it corresponds to caseQ] 



Configuration [2] Type (1, 2) configurations can be generated only by weights vec- 
tors [6, 2] and [4, 4]. They correspond to the cases [2 1 and|2]2. 



Configuration [3] A type (1,3) configuration can be generated only by weights 
vector [4, 2, 2] and it corresponds to case[3] 

Configuration |4] A type (1,4) configuration can be generated only by weights 
vector [2, 2, 2, 2] and it corresponds to case [4] 

Configuration |5] 

Applying Theorem 13.131 configurations of type (2,2) must have one of the 
following weights vectors: 

[5, 1,1,1], [4, 2, 2], [3, 3, 1,1], [2, 2, 2, 2]. 
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• [5, 1, 1, 1] clearly can only generate configuration [5] 1 

• [4, 2, 2] can only results in configurations [5]2 and[5]3. It is necessary to prove 
that they are the only two configurations arising from this weights vector. 

ti i is equal to 4 because it is the max. There are two d.c.'s and for any of 
these d.c.'s the column weight must be even and not zero (if there is a zero 
row then it falls in a smaller configuration). The sum of the column weights 
must be 8 hence, or c-1 and c-2 column weights are 4, or c-1 weights 6 and 
c-2 weights 2 (the case c-1 column weight equal to 2 and c-2 column weight 
equal to 6 is not considered because £x,l = 4). In the first case the only 
possibility is c-l= 4 0, and c-2= 2 2 | , this is configuration 3. In 
the second case c-1 can only be 4 2 , c-2 can be or 2 and that is 
configuration[5]2, or 2 that is the transpose of configuration[5] 3. 

• [3, 3, 1, f] can only results in configurations[5]4 and[3]5. t\ t \ is equal to 3 because it 
is the max. As in the previous case or c-1 and c-2 have weight 4, or c-1 has weight 
6 and c-2 has weight 2. 

i i i 

In the first case, or c- 1 is | 3 I | and c-2 is | 1 3 , or c-1 is 3 1 | and 

c-2 is 3 3 , the first correspond to configuration \5\5 the second to the trans- 
pose of Configuration[5]4. 

In the second case c-1 must be 3 3 , and c-2 must be 1 1 , this is config- 
uration[5]4. 

• [2, 2, 2, 2] evidently can generate only configuration[5]6. 



Configuration [6] 

Applying Theorem 13.131 configurations of type (2,3) must have one of the 
following weights vectors: 

[4, 1, 1, 1, 1] , [3, 3, 1, 1] , [3, 2, 1, 1, 1] , [3, 1, 1, 1, 1, 1] , 
[2, 2, 2, 2], [2, 2, 1,1, 1,1]. 

• [4, 1, 1, 1, 1] is the weights vector for configuration [6)1. This is the only 
possible configurations rising from such weights vector because otherwise 
there would be a column with odd weight. 

• [3, 3, 1, 1] does not result in any valid configurations. In fact all entries are 
odd, hence every d.c. must have two elements to make its column weight 
even, but there are not sufficient entries to fill three columns. 

• [3, 2, 1, 1, 1] can results only in configurations [6]2 and 03. It is necessary 
to prove that they are the only two configurations erasing from this weights 
vector. 

ti i is equal to 3 because it is the max. As discussed previously the column 
weights of each column must be at least 2. The sum of the columns weight 
must be 8 hence there must be one d.c. with weight 4 and two d.c.'s with 
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weight 2. c-1 must have weight 4 since t\ t \ = 3. Considering also that there 
are two d.r.'s and the sum of their row weights must be 8 then it must be 
either that both the row weights are 4 or one is 6 and the other 2. If r-1 has 
weight 6 then t\ t 2 + £l,3 = 3 hence t\ 2 = 2 and £^3 = 1 ( remember that 
ti,2 > h,3), applying Lemma l3"79l c-2 c-3 configuration [6]2 is obtained. If r-1 
has weight 4 then tip + = 1 that implies tip = 1, = (tip > £1,3), 
applying Lemma [3791 configuration |6l 3 is found. 

• [3, 1, 1, 1, 1, 1] cannot generate valid configurations. In fact there are 6 en- 
tries for six positions so the configurations must be full, but this cause both 
d.r.'s to have odd row weight. 

• [2,2,2,2] can only results in configurations [6]4 and[6]6. t\ t \ is equal to 2 
because it is the max. As for the previous case or r-1 has weight 6 and r-2 
has weight 2 or both have weight 4. It is clear that one case generate [6]4 and 
the other [6] 6 or equivalent. 

• [2, 2, 1, 1, 1, 1] can only generate configurations [6]5, in fact it is the only 
possible way to arrange the entries such that Lemma 13.91 is valid for every 
column. 



Configuration [7] 

Applying Theorem 13.131 configurations of type (2, 4) must have one of the 
following weights vectors: 

[2, 2, 1,1, 1,1], [2, 1,1, 1,1, 1,1], [1,1, 1,1, 1,1, 1,1]. 

• [2, 2, 1, 1, 1, 1] can only generate configurations [7J 1 and[7]2. £11 is equal to 
2, and each d.c. must have column weight two because there are four non 
zero columns that must sum to 8. Considering also that there are two d.r.'s 
and the sum of their row weights must be 8 then or one has row weight 6 
and the other 2 or both the row weights are 4. In the first case the d.r. of 
weight 6 must be [2, 2, 1, 1] (the only four entries that sum to six), applying 
Lemma [3791 on all the d.c.'s configuration [7] 1 or equivalent is obtained. In 
the second case both rows must contain the entries 2, 1, 1, (in the proper 
order), in fact if a d.r. contains [2, 2] then the other elements of that row 
must be zeros (otherwise the row weight would be greater than 4) but also 
the elements of the two columns must be zero (column weight must be two). 
This situation cannot be for Lemma l3.101 Once fixed the first row the second 
row can be obtained applying Lemma [3791 on all the d.c.'s and configuration 
[7] 1 or equivalent is obtained. 

• [2, 1, 1, 1, 1, 1, 1] cannot generate valid configurations. In fact £11 = 2 and 
Lemma [3791 c-1 implies that £2,1 = 0. The other six positions must be filled 
with 1-elements but this cause both d.r.'s to have odd row weight. 
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• [1, 1, 1, 1, 1, 1, 1, 1] clearly can only results in configuration[7]3. 
Configuration |8] 

Applying Theorem 13.131 configurations of type (3, 3) must have one of the 
following weights vectors: 

[4, 1, 1, 1, 1] , [3, 2, 1, 1, 1] , [3, 1, 1, 1, 1, 1] , [2, 2, 1, 1, 1, 1] , [2, 1, 1, 1, 1, 1, 1] , [1, 1, 1, 1, 1, 1, 1, 1] 

• [4, 1, 1, 1, 1] cannot generate valid configurations. In fact there are three 
d.c.'s and for any of this column the weight must be even and not zero, the 
only possibility to have a total sum of eight and column weight even is to 
have one column with weight 4 and two with weight 2; the same is true for 
the d.r.'s. This forces the 4-element to be isolated and this cannot be. 

• [3, 2, 1, 1, 1] cannot generate valid configurations. As proved in the previous 
case one d.c. and one d.r. have column weight four and the others weight 
two. is equal to 3, hence the row weight of r-1 and column weight of 
c- 1 must be 4, this implies that the d.m. with the 2-element is isolated, since 

3 10 

it can only be in the row and column with weights 2: 



cannot be. 



1 1 
2 



and this 



of the type 



is obtained. 



[3, 1, 1, 1, 1, 1] can only generate configuration[8]l. As proved in the previous case 
one d.c. and one d.r. have column weight and row weight four and the others weight 
two, moreover = 3. This implies that in r-1 and c-1 there must be only another 
1-element to complete r-1 and c-1 to row and column weight four. Without loss of 
generality it is supposed them to be in t\ i and £2,1- In r-2 and c-2 now there must be 
only another 1-element (column and row weights of the column/rows different from 
the first must have weight 2). If the 1-element falls in £2,2 then there a configuration 
3 1 
1 1 

x 

This configuration is not possible for Lemma [3.10l 
If the 1-element does not fall in £2,2 then there must be a 1-element in position £3^ 
to complete c-2 and one in £2,3 to complete r-2, hence £3 3 = 1 Lemma [3~9l r-3 c-3. 
Configuration^ 1 is obtained. 

[2, 2, 1, 1, 1, 1] can generate only configurations [8] 2 and[8]4. c-1 has weight four, 
since ti t i = 2 and it is not possible to obtain row weight four without using a 2- 
element. Note that it is not possible to fix also r-1 to have weight four because 
the d.r. with weight four could be composed with the other 2-element. Two cases 
are considered. In the first case c-1 is completed with a 2-element, without loss of 
generality it is possible to suppose to be £2.1 = 2. In this case the first column would 
be: \2 2 0| T . In the second case c-1 is completed with two 1-elements, hence the 
first column would be: \2 1 1 | T . 

- c-1 is 1 2 2 0| T . One of r-1, r-2 must have row weight four and it must be com- 



pleted with two 1-elements. Supposing this to be r-1, then for Lemma [3T9| c-2 
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c-3 and the fact that r-2 has already weight 2, £3 2 = £3,3 = 1. Configuration 
[8]2 is obtained. It the d.r. with row weight four is r-2 a configuration that is a 
row permutation of this is obtained. 

- c-1 is |2 1 1| T . If t% t 2 — 2 and ^3 = 0, or tip — and £x,3 = 2 the transpose 
of the previous case is obtained and the resulting configuration is the transpose 
of the previous configuration. Due to Lemma [3~9l r-1 and the fact that the row 
sum cannot be > 4 there are only two possibility for the remaining elements of 
r-1: £1,2 = t\ 3 = 1 or tip = £1,3 = 0. The first case cannot be because there 
is no possibility to place the 2-element satisfying Lemma [3^9l r-2 r-3 c-2 c-3. 
In the second case for Lemma l3.9l applied to r-2 (or r-3) it must be £2(31,2 = 1> 
£2(3), 3 = 2 or £2(3), 2 = 1> £2(3). 3 = 0. It is evident how choosing between 
one of these two possibilities fixes (Lemma [3.9| c-2 c-3) the values in r-3 in a 
configuration that is (apart for row or column permutation) configuration|8]4. 

• [2,1,1,1,1,1,1] can generate only configuration [8] 3. c-1 has weight four, since 
£1.1 = 2 and it is the only 2-element. For the same reason r-1 has row weight 
four, hence tip — £1,3 = £2.1 = £3,1 = 1- Applying Lemma l3~9l r-2. r-3, c-2, c-3 
a configuration that is (apart for row or column permutation) Configuration [8] 3 is 
obtained. 

• [1,1,1,1,1,1,1,1] cannot generate valid configurations. In fact it would be impos- 
sible to have column weight four. 

Configuration [9] 

Applying Theorem 13.131 configurations of type (3, 4) must have one of the 
following weights vectors: 

[2, 2, 1,1, 1,1], [2, 1,1, 1,1, 1,1], [1,1, 1,1, 1,1, 1,1]. 

• [2, 2, 1, 1, 1, 1] cannot generate valid configurations. In fact there are four 
d.c.'s and for any of this column the weight must be at least 2. This implies 
that for each column the weight is exactly 2 otherwise the total sum would 
exceed 8. Also, there are three rows hence, as discussed previously for a 
similar case, the only possibility is to have one row with row weight 4 and 
the other two with row weights 2. If a 2-element lies on a d.r. with row 
weight two that row is complete ans so it is column, hence the 2-element 
is isolated and this is not possible (Lemma [3.101) . It must be that both the 
2-elements lie in the row with row weight four, nothing else lies in the such 
row, but this cannot be because they would be isolated since all d.c.'s have 
column weight two. 

• [2, 1, 1, 1, 1, 1, 1] can only generate our configuration [9jl. t\ t \ is equal to 
2, hence r-1 must have row weight four, otherwise T\ \ would be isolated. 
Hence there are two 1 -elements in r-1. Applying Lemma [3T9l r-2. r-3, c-2, c-3 
a configuration that is (apart for row or column permutation) Configuration 
[9] 1 is obtained. 
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• [1,1,1,1,1,1,1,1] can generate only configuration [9l 2. The row with weight 
four must have all 1-elements, suppose this is the first. Considering Lemma 
13.91 applied to all columns and rows and the fact that row weight of r-2 and 
r-3 are equal to 2, configuration [9]2 or equivalent is obtained. 

Configuration [10] 

Applying Theorem 13.131 configurations of type (4, 4) must have one of the 
following weights vectors: 

[2, 1,1, 1,1, 1,1], [1,1, 1,1, 1,1, 1,1]. 

• [2, 1, 1, 1, 1, 1, 1] cannot generate valid configurations. In fact there are four 
d.c.'s and four d.r.'s, for any of this column the weight must be at least 2. This 
implies that for each column and row the weight is exactly 2 otherwise the 
total sum will exceed 8. Wherever the 2-element lies that row and column is 
completed (every row and column has weight two), hence it is isolated, but 
this cannot be for Lemma l3.10l 

• [1,1,1,1,1,1,1,1] can only generate configuration [9) It is only a case by 
case analysis to prove that any possible distribution of ones, such that form a 
cycle, can be re-conduced to configuration [10] by a combination of rows and 
column permutations. . 

□ 

Proof of theorem 14.101 

Theorem 10.2 (14.101) . Let be M S C m a / 3 7 . The configurations in M that may 
contain a cycles of length 8, are the following^ 



1. 



|C-8|, 
|C-6 C-2|, 
\C-4 C-4|, 
|C-4 C-2 C-2\. 
|C-2 C-2 C-2 C - 21 



For brevity the transposes are omitted 
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C-5 A-l 
A-l A-l 



C-4 A-2 
C-2 



C-4 C-2 
C-2 



C-3 A-l 
C-3 A-l 



C-3 A-l 
A-l C-3 ' 

A-2 A-2 
A - 2 A - 2 ' 

C-4 A-l A-l 
O A-l A-l 



C-3 C-2 A-l 
A-l O A-l 



C-3 A-l O 
A-l A-l C-2 



A-2 C-2 C-2 
C-2 O O 

A-2 A-l A-l 
A-2 A-l A-l 

C-2 A-2 O 
O A-2 C-2 

-2 C-2 A-l A- 
O O A-l A- 
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19. 



20. 



21. 



22. 



23. 



24. 



25. 



26. 



27. 



C-2 A-l A-l O 

O A-l A-l C-2 

A-l A-l A-l A-l 

A-l A-l A-l A-l 



C-3 A-l O 
A-l O A-l 
O A-l A-l 



A - 2 A-l A-l 



C-2 O 



O 



O A-l A-l 



A - 2 A-l A-l 
A-l A-l O 
A-l O A-l 



C-2 O 
A-l A-l 



O 
O 



A-l A-l C-2 



C-2 A-l A-l O 

O A-l O A-l 

O O A-l A-l 

A-l A-l A-l A-l 

A - 1 A - 1 O 

O O A-l A-l 

A-l A-l O O 

A-l O A-l O 

O A-l O A-l 

O O A-l A-l 



Proof. All the configurations appeared in Theorem l3.16l are consider and the rela- 
tive cycle configurations are provided. 
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Configuration Q] gives 

|C — 8] . 

Cycle configuration | J — 8| may be discarded (Lemma Re - 
configuration [2] 

1. Configuration [2l gives 

| C-6 (7-2 | . 

2. Configuration [22 gives 

| C-4 (7-4 | . 

In fact other cycle configurations \C — 6 J — 2|, \J — 6 J — 2|,|J — 6 C — 2| 
and |(7 — 4 J — 4|, | J — 4 J — 4| , | J — 4 (7 — 4| may be discarded because 
in J — 2 and in J — 4 there is no cycle columns (Lemma |4~4U ). 

Configuration [3] gives 

| (7-6 (7-2 (7-2 | . 

Other cycle configurations may be discarded because they contain a d.s. J — 2 or 
J — 6 and in them there is no cycle columns (Lemma |4~4U ). 

Configuration [4] gives 

| (7-2 (7-2 (7-2 (7-2 | . 

Other cycle configurations may be discarded because they contain a d.s. J — 2 and 
in it there is no cycle column (Lemma l4!4l -1). 

Configuration H] 

1. Configuration [5ll gives 

(7-5 A-l 
A - 1 A - 1 " 

Other cycle configurations may be discarded because they contain a d.s. J— 5 
(Lemma 1431). 
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2. Configuration [5] 2 gives 

C-4 C-2 
C-2 

Other cycle configurations may be discarded because they contain a d.s. J— 2 
as the only non-zero matrix in a d.r./d.c. or a J— 4 and that cannot be (Lemma 
O). 

3. Configuration [5j3 gives 

C-4 A-2 

C-2 ' 

Other cycle configurations may be discarded because they contain a d.s. J— 2 
or J — 4 as the only non-zero matrix in a d.r. or d.c. 

4. Configuration [5j4 gives 

C-3 A-l 
C-3 A-l " 

Other cycle configurations may be discarded because they contain a d.s. J— 3 
(Lemma 1431). 

5. Configuration [5j5 gives 

C-3 A-l 
A-l C-3 " 

Other cycle configurations may be discarded because they contain a d.s. J— 3 
(Lemma 1431). 

6. Configuration [5j6 gives 

A-2 A-2 
A-2 A-2 ' 

This cover all the possibilities cycle configuration for the particular case. 
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Configuration [6] 



1. Configuration [6ll gives 



C-4 A 
A 



1 A 
1 A 



Other cycle configurations may be discarded because they contain a d.s. J— 4 
(Lemma [ 



2. Configuration 02 gives 



C 
A 



3 C-2 A 
1 A 



Other cycle configurations may be discarded because they contain a d.s. J— 2 
as the only non-zero matrix in a d.c. or a J — 3 and this cannot be (Lemma 

HjD>. 

3. Configuration [6] 3 gives 

C-3 A- 1 
A-l A-l C-2 

Other cycle configurations may be discarded because they contain a d.s. J— 2 
as the only non-zero matrix in a d.c. or they contain a J — 3 this cannot be 
(Lemma [■ 



A-2 C-2 C-2 
C-2 



4. Configuration [6] 4 gives 



Other cycle configurations may be discarded because they contain a d.s. J— 2 
as the only non-zero matrix in a d.c. or d.r. . 



5. Configuration [6] 5 gives 



A 
A 



2 A 
2 A 



1 A 
1 A 



This cover all the possibilities cycle configuration for the particular case. 
6. Configuration [6] 6 gives 

A-2 C-2 
A-2 C-2 

Other cycle configurations may be discarded because they contain a d.s. J— 2 
as the only non-zero matrix in a d.c. . 
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Configuration |7] 

1. Configuration |7Jl gives 

C-2 C-2 A-l A-l 
A-l A-l 

Other cycle configurations may be discarded because they contain a d.s. J— 2 
as the only non-zero matrix in a d.c. . 

2. Configuration [7]2 gives 

C-2 A-l A-l 
A-l A-l C-2 ' 

Other cycle configurations may be discarded because they contain a d.s. J— 2 
as the only non-zero matrix in a d.c. . 

3. Configuration |7J3 gives 

A-l A-l A-l A-l 
A-l A-l A-l A-l ' 

This cover all the possibilities cycle configuration for the particular case. 



Configuration |8] 

1. Configuration [8jl gives 

C-3 A-l 
A-l A-l . 
A-l A-l 

Other cycle configurations may be discarded because they contain a d.s. J— 3 
(Lemma 1431) . 

2. Configuration [8] 2 gives 

A - 2 A-l A-l 
C-2 

A-l A-l 

Other cycle configurations may be discarded because they contain a d.s. J— 2 
as the only non-zero matrix in a d.r. . 
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3. Configuration [8] 3 gives 



A - 2 A - 1 A - 1 
A - 1 A - 1 
A - 1 A - 1 

This cover all the possibilities cycle configuration for the particular case. 

4. Configuration [8]4 gives 

C-2 
A - 1 A - 1 
A-l A-l C-2 

Other cycle configurations may be discarded because they contain a d.s. J— 2 
as the only non-zero matrix in a d.r. or d.c. . 



Configuration [9] 

1. Configuration [9)1 gives 

A-l A-l A-l A-l 
A - 1 A - 1 
A - 1 A - 1 

This cover all the possibilities cycle configuration for the particular case. 

2. Configuration [9] 2 gives 

C-2 A-l A-l 
A-l A - 1 . 
A - 1 A - 1 

Other cycle configurations may be discarded because they contain a d.s. J— 2 
as the only non-zero matrix in a d.c. . 



Configuration [10] 

A - 1 A - 1 

A-l A-l 

A-l A - 1 ' 

A - 1 A - 1 

This cover all the possibilities cycle configuration for the particular case. 

For Remark 13.71 it is not necessary to study the cycle configurations that are 
transposed of the one considered. Hence it has been proved that the listed configu- 
rations are the only valid. □ 
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Proof of theorem 15.161 

Theorem 10.3 (15.161) . Let be M € C miQ , )( g i7 . The configurations in M that may 
contain a cycles of length exactly 8, are the following^ 

1. 

|C-8|, a(p) = m/4 

2. 

&-h J 2 -I 
J 3 -l A 4 -l ' 

e(p 1 )-e(p 2 )-e(p 3 )+e(p*) = ±2s(p 1 ), 

3. 

C 1 - 4 J 2 - 2 
C 3 - 2 ' 

±s(p 3 ) = 2s(p 1 ) 

4. 

C 1 -?, J 2 -I 
J 3 -l C 4 -3 ' 

e (pl) _ e (p2) _ e (p3) + e(p 4) = ±a(p l) ± s(p 3 }) 

5. 

A 1 - 2 A 2 - 2 
A 3 -2 A 4 -2 ' 

e(p l) + e(p l) _ e (p2) _ e(p 2) _ e(p 3) _ £(p 3 } + e(p4) + e(p4) = ^ 

6. 

C 1 - 4 J 2 - 1 J 3 - 1 
O A 4 -l A 5 -l ' 

e(p 2 )-e(p 3 )-e(p 4 ) + e(p 5 )^±2s(p 1 ), 

7. 

C 1 - 3 O J 2 - 1 
J 3 -l C A -2 J 5 -l ' 

£(p l) _ e (p2) _ e(p 3) + £(p5) = ±s(p l) ± s(p 4) 9 

9 Configurations with two or more weight-2 circulants in the same row or column are n not listed 
since they always contain a cycle of at most 8 ("Lemma [5,15b . hence they do not add any information. 
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A 1 - 2 A 2 - 1 A 3 - 1 
A 4 -2 A 5 -l A 6 -l ' 

e(p l) + e(p l) _ e(p 2) _ e(p 3) _ e( /) _ e( /) + e(p 5) + e(p 6) = , or 
e(p 2) _ e(p 3) + e(p 4) _ e(p 4) _ e(p 5 } + £(p6) = ±s(p l) j or 
e(p 2) _ e( p3) + e(p l) _ e(p l) _ e(p 5) + e( /) = ±s(p 4) . 

C 1 -2 J 2 -2 O 
O J 3 -2 C 4 -2 ' 

C 1 - 2 O J 3 - 1 J 4 - 1 
O C 2 -2 J 5 -l J 6 -l ' 

e(p 3) _ e( p4) _ e(p 5) + e(p 6) = ±s(p l) ± s( p2 }) 

A 4 -l A 2 -l A 3 -l A 4 -l 
A 5 -l A 6 -l A 7 -l A 8 -l ' 

e(p x ) + e(p 2 ) - e(p 3 ) - e(p 4 ) - e(p 5 ) - e(p 6 ) + e(p 7 ) + e(p 8 ) = , or 
tip 1 ) - e(p 2 ) + e(p 3 ) - e(p 4 ) - e(p 5 ) + e(p 6 ) - e{p 7 ) + e(p 8 ) = , or 
e(p l) _ e(p 2 ) _ e(p 3 ) + e (p4) _ e( p5) + e( /) + _ e(p 8) = o . 



C 1 - 3 J 2 - 1 O 
J 3 -l O A 4 -l , 
O A 5 -l A 6 -l 

e(p l) _ e(p 2 ) _ e(p 3 ) + e (p4) + e( p5) _ £(p 6) = ±s(p l )> 



ji _ 2 A 2 - 1 A 3 - 1 
C 4 -2 O O 
O A 5 -l A 6 -l 

e(p 2 )-e(p 3 )-e(p 5 ) + e(p 6 ) = ±s(p 4 ) 
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14. 



A 1 - 2 A 2 - 1 A 3 - 1 
A 4 -l A 5 -l O 
A 6 - 1 O A 7 - 1 

tip 1 ) + tip 1 ) ~ e(p 2 ) - e(p 3 ) - e(p 4 ) + e(p 5 ) - e(p 6 ) + e(p 7 ) = , or 

e(p 2 ) ~ e(p 3 ) + e(p 4 ) - e(p 5 ) - e(p 6 ) + e(p 7 ) = i^p 1 ), 



75. 



C 1 - 2 
J 2 
J 4 



o o 

1 A 3 - 1 O 
1 J 5 - 1 C 6 - 2 



id 



77. 



18. 





A 1 - 


2 A 2 -l A 3 -l 


A 4 -l 






A 5 - 


1 A 6 -l O 





1 







O A 7 -l 


A 8 -l 




e (p2) _ e (p3) _ 


e(p 4 )+e(p 5 )+e(p 6 


) " <P T ) - 






c 1 - 


2 J 2 - 1 J 3 - 1 












A 4 - 1 O 


A 5 - 1 









O A 6 - 1 


A 7 -l 




e(p 2) _ e (p3) _ 


e(p 4 ) + e(p 5 ) + e(p 6 


) " <P 7 ) - 






A 1 - 


1 A 2 -l O 


O 






A 3 - 


1 A 4 - 1 












A 5 - 1 O 


A 6 - 1 


1 







O A 7 -l 


A 8 -l 




eCp 1 ) - e(p 2 ) 


- <P Z 


) + e(p 4 ) + e(p 5 )-6 


(p 6 )-e(p 7 )+e(p 8 



0. 



Proof. As for the theorems 15.31 and 15.71 the proof is once more divided in smaller 
lemmas each considering a particular case. 

An 8-cycle is formed by four cycle comumns called y,z,v,u and four cycle rows 
called x,t,w,l. In every 8-cycle there are 8 cycle points, they take the name of the 
column and row where they lie: (x,y), (x,z), (t,y), (t,v), (w,u), (w,z), (l,u), (l,v). 
The notation is used to compute the conditions. 
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Lemma 10.4. There is a 8-cycle and there are no 6-cycles or 4-cycles in case\l\if 
and only if 

A\mand s(p)=m/4. 

Proof. There is a 8-cycle and no 6-cycle or 4-cycle if and only if g = 6. Applying 
Prop. 15-11 to the case g = 8 and M = C, this is equivalent to 

m 

8 = g= " 



cd(m, s) 



i.e. m = 4gcd(m,s). In particular, 4|m and m/4|s, but s < m/2, so that s = 
m/4. On the other hand, if m is divisible by 4 and s = m/4 then gcd(m, s) = s = 
m/4. □ 

Lemma 10.5. There is a 8-cycle in case\I\if and only if 

e(p 1 )-e(p 2 )-e(p 3 )+e(p 4 ) = ±2 S (p 1 ) 

Proof. It can be assumed that there is a 8-cycle if and only if, simultaneously, 
cycle columns y and z lie in C , cycle point (t,v) lies in C , cycle point (w,u) 
lies in A 2 , cycle point (I, v) lies in A 3 and cycle point (I, u) lies in A 4 . Applying 
Proposition 12.91 -4 to cycle column y and cycle column z yields 

(46) x - t = ±s 1 , 

(47) x-w = =Fs x . 

Since cycle point (t, v) lies in C , applying Proposition 12.91 1 

(48) v = t + e^ 1 ) . 

Since cycle point (w, u) lies in A 2 , applying Proposition 12.91 -1 

(49) u = w + e(p 2 ) . 
Since cycle point (I, v) lies in A 3 , applying Lemma [3^41 

(50) v = l + e(p 3 ) . 
Since cycle point (I, u) lies in A 4 , applying Lemma [3~4l 

(51) u = l + e{p 4 ). 
From g8]l, dH, and g5>, (ED : 

(52) / + e(p 3 ) ^t + e^ 1 ). 

(53) / + e(p 4 ) =w + e(p 2 ). 
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Subtracting d52J) from (53) : 

(54) e{p 4 ) - e(p 3 ) = w - t + e{p 2 ) - e{p l ) 

(55) e(p 2 ) - e{p l ) - e(p 3 ) + e(p 4 ) = w-t 

but (to — i) = (x — t) — (x — w) hence from (l46l ) and (|47"1) . 

(56) w-t = ±2s(p 1 ) 

The condition can eb obtained from (l54l and (l56l >. □ 

Lemma 10.6. There is an 8-cycle in case\3\if and only if 

e(p 2 )-e(p 2 ) = ±s(p 3 )±2s(p 1 ) 

Proof. It can be assumed that there is a 8-cycle if and only if columns y and z lie in 
C 1 , cycle row / lies in C 3 , cycle point (t, v) lies in A 2 , and cycle point (u>, u) lies 
in A 2 . Applying Proposition 12.91 4 to cycle column y and cycle column z yields 

(57) x - t = ±s 1 , 

(58) x - w = ^s 1 . 
Since cycle point (t, v) lies in A 2 , applying Lemma l3~4l 

(59) v = t + e(p 2 ) . 

Since cycle point (w, u) lies in A 2 , applying Lemma [341 

(60) u = w + e(p 2 ) . 
Since cycle row I lies in C 3 , applying Proposition 12.91 3 

(61) v-u = ±s 3 . 

Subtracting ((60]> to ((59]) and substituting the result into (J6B and using ((57]>, (l58l) 
the desired result is obtained. □ 

Lemma 10.7. There is a 8-cycle in case^if and only if 

^P 1 ) ~ <P 2 ) ~ <p 3 ) + e(p 4 ) = ±s(p 1 ) ± s(p 4 ) 

Proof. It can be assumed that there is a 8-cycle if and only if cycle row x lies in 
C 1 , cycle row I lies in C 4 , cycle point (w, z) lies in C 1 , cycle point (w, u) lies in 
J 2 , cycle point (i, y) lies in J 3 , and cycle point (t, v) lies in C 4 . Since cycle row 
x lies in C 1 , applying Proposition 12.91 -3 

(62) z - y = ±s 1 . 
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Since cycle row I lies in C 4 , applying Proposition 12. 91 3 

(63) u - v = ±s 4 . 

Since cycle point (w, z) lies in C , applying Lemma [3~4l 

(64) z^w + eip 1 ). 
Since cycle point (w, u) lies in J 2 , applying Lemma l3~4l 

(65) u = w + e(p 2 ) . 
Since cycle point (£, y) lies in J 3 , applying Lemma l3~4l 

(66) y = t + e(p 3 ) . 
Since cycle point (t, v) lies in C 4 , applying Lemma [3~4l 

(67) v = t + e(p 4 ) . 

The desired result can be computed form the previous formulas. □ 
Lemma 10.8. There is a 8-cycle in case\5}if and only if 

e (pl) + e (pl) _ e (p2) _ e (p2) _ e (p3) _ e (p3) + e(p 4) + e(p 4) = , 

Proof. Two cases are considered, the first when there is a circulant column or row 
in A 1 , the other when there is not. The aim is now to prove that in the first case 
no 8 cycle can be formed for the class of quasi-cyclic codes considered. Thanks 
to the equivalence of the configurations over row column permutation it is possible 
to suppose that the cycle column/row lies in A 1 , moreover only the case of cycle 
column need to be considered since the case with A 1 containing a cycle row is the 
transpose of this. 

To have a cycle column in A 1 it must be a weight 2 circulant. For lemma l5~15l 
A 2 and A 3 cannot be weight-2 circulants. The two cycle points lying in A 3 cannot 
form a cycle column because A 3 is a weight- 1 circulant, hence they must be part of 
two different cycle columns but this is not possible because in C 1 there are no other 
points to form cycle columns. This shows how cycles exist in this configuration if 
and only if none of the four circulants contain a cycle column or a cycle row. 

If there are no cycle columns or row in any circulant it can be assumed that 
there is a 8-cycle if and only if cycle point (x, y) lies in A 1 , cycle point (I, u) lies 
in A 1 , cycle point (x, z) lies in A 2 , cycle point (I, v) lies in A 2 , cycle point (w, u) 
lies in A 3 , cycle point (t, y) lies in A 3 , cycle point (w, z) lies in A 4 , cycle point 
(t, v) lies in A 4 , since cycle point (x, y) lies in A 1 , applying Lemma [3~4l 

(68) y = x + e(/) . 
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Since cycle point (I, u) lies in A , applying Lemma [3~4l 

(69) u = l + e(p l ). 
Since cycle point (x,z) lies in A 2 , applying Lemma l3~4l 

(70) z = x + e(p 2 ) . 
Since cycle point (I, v) lies in A 2 , applying Lemma l3~4l 

(71) v = l + e{p 2 ). 

Since cycle point (w, u) lies in A 3 , applying Lemma l3~4l 

(72) u = w + e(p 3 ) . 
Since cycle point (t, y) lies in A 3 , applying Lemma l3~4l 

(73) y = t + e(p 3 ) . 

Since cycle point (w, z) lies in A 4 , applying Lemma l3~4l 

(74) z = w + e(p A ) . 
Since cycle point (t, v) lies in A 4 , applying Lemma l3~4l 

(75) v = t + e{p A ) . 

From (|69l and d72l eq.|76]is obtained, substituting (in order) (1711 ). (T75T >. (1731 . 
d70l and d74l the desired condition is found. 



(76) l-w + e{p l ) - e(p 3 ) = 

Remark 10.9. x In the case when all the circulants in the configuration are weight- 
1 circulants then this condition is redundant. In fact is such case the condition 
became: 

2e{p 1 ) - 2e{p 2 ) - 2e{p 3 ) + 2e(p 4 ) = 

that is equivalent to case\3\ofTheorem \5.3\ hence this condition assure the existence 
of a A-cycle and not of an 8-cycle. 

□ 

Lemma 10.10. There is a 8-cycle in case^if and only if 

e(p 2 )-e(p 3 )-e(p 4 ) + e(p 5 ) = ±2 S (p 1 ) 
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Proof. It can be assumed that there is a 8-cycle if and only if cycle columns y and 
z lie in C 1 , cycle point (t, v) lies in A 2 , cycle point (to, u) lies in A 3 , cycle point 
(l,v) lies in A 4 , and cycle point (l,u) lies in A 5 . Applying Proposition 12.91 4 to 
cycle column y and cycle column z yields 

(77) x - t = ±s 1 , 

(78) x - w = =fs x . 
Since cycle point (t, v) lies in A 2 , applying Lemma l3~4l 

(79) v = t + e{p 2 ) . 

Since cycle point (w, u) lies in A 3 , applying Lemma [3~4l 

(80) u = w + e{p 3 ) . 
Since cycle point (I, v) lies in A 4 , applying Lemma |3~4| 

(81) v = l + e{p A ). 
Since cycle point (I, u) lies in A 5 , applying Lemma [3^41 

(82) u = l + e(p 5 ) . 
Substituting into G3) : 

(83) w-t = ±2s l . 

Substituting d79]), ([80]>, ([SB and ((82]> in d83j the desired condition is obtained. □ 
Lemma 10.11. There is a 8-cycle in case\7\if and only if 

e (pl) _ e (/) _ e(p 3 ) + e (pB) = ±s(p l) ± 

Proof. It can be assumed that there is a 8-cycle if and only if cycle columns y lies 
in C , cycle columns u lies in C 4 , cycle point (t, u) lies in C , cycle point (x, z) 
lies in J 2 , cycle point (7, lies in J 3 , and cycle point (w, z) lies in J 5 . Applying 
Proposition 12.91 -2 to cycle column y 

(84) x - t = ±s 1 . 
Applying Proposition 12.91 -2 to cycle column z 

(85) w - I = ±s 2 . 
Since cycle point (t, v) lies in C 1 , applying Lemma [3^41 

(86) v = t + eOp 1 ) . 
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Since cycle point (x, z) lies in J 2 , applying Lemma [3~4l 

(87) z = x + e(p 2 ) . 
Since cycle point (l,v) lies in J 3 , applying Lemma [3~4l 

(88) v = I + e{p 3 ) . 
Since cycle point (w, z) lies in J 5 , applying Lemma l3~4l 

(89) z = w + e(p 5 ) . 

Substituting in the proper order the desired result is obtained. □ 

Lemma 10.12. There is a 8-cycle in case\8\if and only if 

e(p l ) + e^ 1 ) - e(p 2 ) - e(p 3 ) - e(p 4 ) - e(p A ) + e(p 5 ) + e(p 6 ) = , or 
e(p 2 ) - e(p 3 ) + e(p A ) - e(p A ) - e(p 5 ) + e(p 6 ) = ±s(p 1 ) , or 
e(p 2 ) - e(p 3 ) + eip 1 ) - e(p l ) - e(p 5 ) + e(p 6 ) = ±s(p 3 ) . 

Proof. Two cases must be considered, the first when nor A 1 neither A 4 contain a 
cycle column or a cycle row, the second when one of them contains a cycle row. 
Note how for lemma [5TT51 A 1 and A 4 cannot be weight-2 circulants at the same 
time hence they cannot both contain a cycle row or cycle column. The situation 
with only one cycle columns in A 1 or in A 4 can be dismissed, in fact it can be 
assumed the cycle column to lie in A 1 . Then the two cycle points in A 4 must be in 
two different columns but these columns cannot exist because there are not other 
point in A 1 to complete them. 

• If there are no cycle columns or row in any circulant it can be assumed that 
there is a 8-cycle if and only if cycle point (x, y) lies in A 1 , cycle point (I, u) 
lies in A 1 , cycle point (x, z) lies in A 3 , cycle point (I, v) lies in A 2 , cycle 
point (w, u) lies in A 4 , cycle point (t, y) lies in A 4 , cycle point (w, z) lies 
in A 6 , and cycle point (i, v) lies in A 5 . Since cycle point (x, y) lies in A 1 , 
applying Lemma [3~4l 



(90) y = x + e(p l ). 

Since cycle point (l,u) lies in A 1 , applying Lemma [3~4l 

(91) u = I + eip 1 ) . 

Since cycle point (x, z) lies in A 3 , applying Lemma [3~4l 

(92) z = x + e(p 3 ) . 

Since cycle point (I, v) lies in A 2 , applying Lemma [341 

(93) v = I + e(p 2 ) . 
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Since cycle point (w, u) lies in A 4 , applying Lemma [3~4l 

(94) u = w + e(p 4 ) . 

Since cycle point (t, y) lies in A 4 , applying Lemma [3^41 

(95) y = t + e(p 4 ) . 

Since cycle point (w, z) lies in A 6 , applying Lemma l3~4l 

(96) z = w + e(p 6 ) . 

Since cycle point (t, v) lies in A 5 , applying Lemma l3~4l 

(97) v = t + e(p 5 ) . 

From d£B and d94j) eq. 198 1 is obtained and substituting (in order) (1931 ), d97| ), 
|), d90]>, ([92j> and ([96l) the listed condition is obtained. 



(98) I - w + e(p x ) - e(p 4 ) = 

• If there is a cycle row in a circulant it can be assumed that it is in C 1 . Then 
there is not any other cycle row or cycle column in C 4 , hence it can be 
assumed that there is a 8-cycle if and only if cycle row x lies in C 1 , cycle 
point (w, z) lies in A 4 , cycle point (t, y) lies in A 4 , cycle point (t, v) lies 
in A 6 , cycle point (w, u) lies in A 5 , cycle point (I, u) lies in A 2 , and cycle 
point (l,v) lies in A 3 . Since cycle row x lies in C 1 , applying Proposition 
1251 3 

(99) y - z = ±s 1 . 

Since cycle point (w, z) lies in A 4 , applying Lemma [3^41 

(100) z = w + e(p 4 ) . 

Since cycle point (t, y) lies in A 4 , applying Lemma 1341 

(101) y = t + e{p 4 ). 

Since cycle point (t, v) lies in A 6 , applying Lemma [3^41 

(102) v = t + e(p 6 ). 

Since cycle point (w, u) lies in A 5 , applying Lemma [3~4l 

(103) u = w + e(p 5 ). 

Since cycle point (/, u) lies in A 2 , applying Lemma l3~4l 

(104) u = I + e(p 2 ) . 
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Since cycle point (I, v) lies in A 3 , applying Lemma [3~4l 
(105) v = I + e(p 3 ) . 

Staring from (g9]l and substituting in order (flOOl . (fTOTT) . (fT0"2l . (fT03l> . (fT04l) 
and (11051 ) the result is obtained: 

e(p 2 ) ~ e(p 3 ) + e(p 4 ) - e(p 4 ) - e(p 5 ) + e(p 6 ) = isfr 1 ) 

The last condition 

e(p 2 ) - e(p 3 ) + eip 1 ) - eip 1 ) - e(p 5 ) + e(p 6 ) = ± S (p 3 ) 
Is equivalent to this case but with cycle row laying in C 4 instead of C 1 . 

□ 

Lemma 10.13. There is a 8-cycle in case^if and only if 

e(p 2 ) ~ e(p 2 ) + e(p 3 ) - e(p 3 ) = ±s(p 1 ) ± S (p 4 )- 

Proof. It can be assumed that there is a 8-cycle if and only if cycle column y lies 
in C 1 , cycle column u lies in C 4 , cycle point (x, z) lies in A 2 , cycle point (t, v) 
lies in A 2 , cycle point (w, z) lies in A 3 , and cycle point (I, v) lies in A 3 . Applying 
Proposition 12.91 -2 to cycle column y 

(106) x - t = ±s 1 . 
Applying Proposition 12.91 -2 to cycle column u 

(107) w - I = ±s 4 . 

Since cycle point (x, z) lies in A 2 , applying Lemma [3~4l 

(108) z = x + e(p 2 ) . 
Since cycle point (t, v) lies in A 2 , applying Lemma l3~4l 

(109) v = t + e(p 2 ) . 

Since cycle point (w, z) lies in A 3 , applying Lemma [3~4l 

(110) z = w + e(p 3 ). 
Since cycle point (I, v) lies in A 3 , applying Lemma [3^41 

(111) v = l + e(p 3 ). 

Substituting (fT08l) and (fT09l into (TT06T> and then (fTOTT) . (fTTOl and (EB, the desired 
result is obtained. □ 
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Lemma 10.14. There is a 8-cycle in case\W\if and only if 

e (p3) _ e (p4) _ e(p 5 } + e(p6) = ±s(p l } ± s(p2) 

Proof. This can be prove with identical reasoning used to prove Lemma [T0.14l □ 
Lemma 10.15. There is a 8-cycle in case\T7}if and only if 

<P l ) + <P 2 ) ~ <p 3 ) - e(p 4 ) - e(p s ) - e(p 6 ) + e(p 7 ) + e(p s ) = , or 
e(p l ) - e(p 2 ) + e(p 3 ) - e(p 4 ) - e(p 5 ) + e(p 6 ) - e(p 7 ) + e(p 8 ) = , or 
tip 1 ) - e(p 2 ) - e(p 3 ) + e(p 4 ) - e(p 5 ) + e(p 6 ) + e(p 7 ) - e(p 8 ) = . 

Proof. To prove this lemma it is sufficient to consider that there is one cycle point 
for each circulant in the configuration, hence even if the circulant may have weight 
2 they act as weight one. Such situation is identical to the class of codes used by 
Fossorier in ll28l . Theorem 14. 1 1 presented in section |4] can be applied directly to 
this class. The three listed conditions are the expansion of the theorem apply to the 
only possible cycles in such configuration. 

In particular fixed the first cycle column to be fixed between Ai and A5 there are 
three possible column that the second row can lie in. For each of this possibilities 
there are two possible positions of the remaining two cycle columns. A total of 
six possible cycles are obtained but they are two by two equivalent by row swap 
(seel26l). Since there are only two rows a row swap changes the sign to all the terms 
hence only three conditions are required. 




Figure 26: All possible 8-cycles for configuration 11 

□ 

Lemma 10.16. There is a 8-cycle in case\l2\if and only if 

e (pl) _ e (p2) _ e (p3) + e (p4) + e(p 5) _ e(p 6) = ±a(p l) . 

Proof. It can be assumed that there is a 8-cycle if and only if cycle column y lies in 
C 1 , cycle point (t, v) lies in C 1 , cycle point (re, z) lies in A 2 , cycle point (I, v) lies 
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in A 3 , cycle point (I, u) lies in A 4 , cycle point (w, z) lies in A 5 , and cycle point 
(w, u) lies in A 6 . Applying Proposition 12.91 -2 to cycle column y 

(112) x-t = ±s 1 . 
Since cycle point (t, v) lies in C 1 , applying Lemma [3^41 

(113) v^t + eip 1 ). 

Since cycle point (x, z) lies in A 2 , applying Lemma l3~4l 

(114) z = x + e{p 2 ). 
Since cycle point (I, v) lies in A 3 , applying Lemma [3^41 

(115) v = l + e(p 3 ). 
Since cycle point (I, u) lies in A 4 , applying Lemma [3^41 

(116) lem : u = I + e(p 4 ) . 
Since cycle point (w, z) lies in A 5 , applying Lemma l3~4l 

(117) z = w + e(p 5 ). 

Since cycle point (w, u) lies in A 6 , applying Lemma [341 

(118) u = w + e(p e ). 

Starting from (11121 ) and substituting the other equations it is easy to find the desired 
result. □ 

Lemma 10.17. There is a 8-cycle in case\l3\if and only if 

e(p l) _ e (pl) + e(p 2) _ e(p 3) _ e (p5) + e (p6) = ±s(p 4) . 

Proof. It can be assumed that there is a 8-cycle if and only if cycle row x lies in 
C , cycle point (w, z) lies in A 1 , cycle point (t, y) lies in A 1 , cycle point (£, w) lies 
in A 2 , cycle point (u;, u) lies in A 3 , cycle point (I, v) lies in A 5 , and cycle point 
(I, u) lies in A 6 . Since cycle row x lies in C 4 , applying Proposition 12.91 3 

(119) z-y = ±s 4 . 

Since cycle point (w, z) lies in A 1 , applying Lemma [3~4l 

(120) z = w + e(p l ). 
Since cycle point (t, y) lies in A 1 , applying Lemma [3^41 

(121) y = t + e (pi). 
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Since cycle point (t, v) lies in A 2 , applying Lemma [3^41 

(122) v = t + e(p 2 ) . 

Since cycle point (w, u) lies in A 3 , applying Lemma [3~4l 

(123) u = w + e(p 3 ) . 
Since cycle point (I, v) lies in A 5 , applying Lemma l3~4l 

(124) v = I + e(p 5 ) . 
Since cycle point (I, u) lies in A 6 , applying Lemma l3~4l 

(125) u = I + e(p 6 ) . 

Starting from (11191 ) and substituting the other equations it is easy to find the desired 
result. □ 

Lemma 10.18. There is a 8-cycle in case\14\if and only if 

t{p l ) + eCp 1 ) - e{p 2 ) - e{p 3 ) - e(p 4 ) + e(p 5 ) - e(p 6 ) + e(p 7 ) = , or 

e (p2) _ e (p3) + e(p 4) _ _ £(p 6) + £(p 7) = ±s(p l) 

Proof. Two possible cases must be considered. The first is when the two point 
laying in A 1 form a cycle column (or row) the second when they do not. 

• If there is a cycle column in A 1 it can be assumed that there is a 8-cycle if 
and only if cycle column y lies in A 1 , cycle point (x, z) lies in A 3 , cycle 
point (t, v) lies in A 2 , cycle point (/, v) lies in A 5 , cycle point (I, u) lies in 
A , cycle point (w, z) lies in A , and cycle point (w, u) lies in A 6 . Applying 
Proposition 12.91 2 to cycle column y 



(126) x - t = ±s 1 . 

Since cycle point (x, z) lies in A 3 , applying Lemma [341 

(127) z = x + e(p 3 ). 

Since cycle point (t, v) lies in A 2 , applying Lemma [3^41 

(128) v = t + e{p 2 ). 

Since cycle point (I, v) lies in A 5 , applying Lemma [3~4l 

(129) v = l + e(p 5 ). 

Since cycle point (/, u) lies in A 4 , applying Lemma l3~4l 

(130) u = I + e(p 4 ) . 
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Since cycle point (w, z) lies in A 7 , applying Lemma [3~4l 

(131) z = w + e(p 7 ). 

Since cycle point (w, u) lies in A 6 , applying Lemma l3~4l 

(132) u = w + e(p 6 ). 

Starting from (11261 ) and substituting the other equations it is easy to find the 
condition in the statement. 

If the two points in A 1 form a cycle row and not a cycle column it is sufficient 
to transpose the matrix (it is a symmetric 3x3 matrix) and the same condition 
is valid. 

• If there is no cycle column or row in A 1 it can be assumed that there is a 
8-cycle if and only if cycle point (x, y) lies in A , cycle point (I, u) lies in 
A 1 , cycle point (x,z) lies in A 2 , cycle point (l,v) lies in A 3 , cycle point 
(w, u) lies in A 4 , cycle point (w, z) lies in A 5 , cycle point (t, y) lies in A 6 , 
and cycle point (t, v) lies in A . Since cycle point (x, y) lies in A 1 , applying 
Lemma l3.4l 

(133) y^x + eip 1 ). 

Since cycle point (/, u) lies in A 1 , applying Lemma 1341 

(134) u^l + eip 1 ). 

Since cycle point (x, z) lies in A 2 , applying Lemma l3~4l 

(135) z = x + e(p 2 ). 

Since cycle point (/, v) lies in A 3 , applying Lemma l3~4l 

(136) v = l + e(p 3 ). 

Since cycle point (w, u) lies in A 4 , applying Lemma [341 

(137) u = w + e(p 4 ). 

Since cycle point (w, z) lies in A 5 , applying Lemma l3~4l 

(138) z = w + e(p 5 ). 

Since cycle point (t, y) lies in A 6 , applying Lemma 1341 

(139) y = t + e(p 6 ). 

Since cycle point (t, v) lies in A 7 , applying Lemma l3~4l 

(140) v = t + e{p 7 ). 

Starting from (11341 ) and (11371 ) eq. U41l is obtained and substituting the other 
equations into this the desired condition is found. 

(141) l + e{p l ) = w + e{p A ). 
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□ 

Lemma 10.19. There is a 8-cycle in case\l5\if and only if 

e(p 2) _ e (p3) _ e(p 4) + e (pB) = ±s(p l) ± s(p 6) . 

Proof. It can be assumed that there is a 8-cycle if and only if cycle row x lies in 
C 1 , cycle column v lies in C 6 , cycle point (to, z) lies in A 2 , cycle point (to, u) lies 
in A 3 , cycle point (t, y) lies in A 4 , and cycle point (I, u) lies in A 5 . Since cycle 
row x lies in C 1 , applying Proposition 12.91 3 

(142) z-y = ±s 1 . 
Applying Proposition 12.91 -2 to cycle column v 

(143) t-l = ±s 6 . 

Since cycle point (to, z) lies in A 2 , applying Lemma 1341 

(144) z = w + e(p 2 ) . 

Since cycle point (w, u) lies in A 3 , applying Lemma l3~4l 

(145) u = w + e(p 3 ) . 
Since cycle point (t, y) lies in A 4 , applying Lemma l3~4l 

(146) y = t + e(p 4 ) . 
Since cycle point (I, u) lies in A 5 , applying Lemma l3~4l 

(147) u = l + e(p 5 ). 

Starting from (| 142b and substituting the other equations it is easy to find the desired 
result. □ 

Lemma 10.20. There is a 8-cycle in case\l6\if and only if 

e (p2) _ e(p 3 } _ e(p 4 } + e(p5) + e(p6) _ e(p7) = ±a(p l )? 

Proof. It can be assumed that there is a 8-cycle if and only if cycle point (x, y) lies 
in A 1 , cycle point (t, y) lies in A 5 , cycle point (x, z) lies in A 4 , cycle point (t, v) 
lies in A 6 , cycle point (to, z) lies in A 8 , cycle point (to, u) lies in A 7 , cycle point 
(I, v) lies in A 2 , and cycle point (I, u) lies in A 3 . 

Since cycle point (x, y) lies in A , applying Lemma [3741 

(148) y = x + e(p 1 ) . 
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Since cycle point (t, y) lies in A 5 , applying Lemma [3~4l 

(149) y = t + e(p 5 ) . 
Since cycle point (x, z) lies in A 4 , applying Lemma l3~4l 

(150) z = x + e(p 4 ) . 
Since cycle point (t, v) lies in A 6 , applying Lemma [3T4l 

(151) v = t + e{p 6 ). 

Since cycle point (w, z) lies in A 8 , applying Lemma l3~4l 

(152) z = w + e(p 8 ). 

Since cycle point (w, u) lies in A 7 , applying Lemma [3~4l 

(153) u = w + e(p 7 ) . 
Since cycle point v) lies in A 2 , applying Lemma [3~4l 

(154) v = l + e(p 2 ) . 
Since cycle point u) lies in A 3 , applying Lemma l3~4l 

(155) u = I + e(p 3 ) . 

Starting from (11481 ) and substituting the other equations it is easy to find the desired 
result. □ 

Lemma 10.21. There is a 8-cycle in case\T7\if and only if 

e(p 2 ) ~ e(p 3 ) - e(p A ) + e(p 5 ) + e(p 6 ) - e(p 7 ) = ±s(p l ) . 

Proof. It can be assumed that there is a 8-cycle if and only if cycle column y lies 
in C 1 , cycle point (x, z) lies in A 2 , cycle point (t, v) lies in A 3 , cycle point (w, z) 
lies in A , cycle point (w,u) lies in A 5 , cycle point (l,v) lies in A 6 , and cycle 
point (7, n) lies in A 7 . Applying Proposition 12.91 2 to cycle column y 

(156) x-t = ±s 1 . 

Since cycle point (x, z) lies in A 2 , applying Lemma l3~4l 

(157) z = x + e{p 2 ). 
Since cycle point (t, v) lies in A 3 , applying Lemma [3^41 

(158) v = t + e(p 3 ). 
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Since cycle point (w, z) lies in A 4 , applying Lemma [3~4l 

(159) z = w + e(p 4 ) . 

Since cycle point (w, u) lies in A 5 , applying Lemma 13741 

(160) u = w + e(p 5 ) . 
Since cycle point (I, v) lies in A 6 , applying Lemma 13741 

(161) v = l + e(p 6 ). 
Since cycle point (I, u) lies in A 7 , applying Lemma 13741 

(162) u = I + e(p 7 ) . 

Starting from (11561 ) and substituting the other equations it is easy to find the desired 
result. □ 

Lemma 10.22. There is a 8-cycle in case\18\if and only if 

e (pl) _ e (p2) _ e (p3) + e(p 4) + £(p 5) _ £(p 6) _ + £(p 8) = Q . 

Proof. It can be assumed that there is a 8-cycle if and only if cycle point (x, y) lies 
in A 1 , cycle point (a?, z) lies in A 2 , cycle point (t, y) lies in A 3 , cycle point (t, v) 
lies in A 4 , cycle point (w, z) lies in A 5 , cycle point (w, u) lies in A 6 , cycle point 
(I, v) lies in A 7 , and cycle point (7, u) lies in A 8 . Since cycle point (x, y) lies in 
A 1 , applying Lemma [3741 

(163) y = £ + 

Since cycle point (x,z) lies in A 2 , applying Lemma l3~4l 

(164) z = x + e(p 2 ) . 
Since cycle point (t, y) lies in A 3 , applying Lemma [3~4l 

(165) y = t + e(p 3 ). 
Since cycle point (t, v) lies in A 4 , applying Lemma [3741 

(166) v = t + e(p 4 ) . 

Since cycle point (w, z) lies in A 5 , applying Lemma 13741 

(167) z = w + e(p 5 ). 

Since cycle point (w, u) lies in A 6 , applying Lemma [3741 

(168) u = w + e(p 6 ). 
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Since cycle point (I, v) lies in A 7 , applying Lemma [3741 

(169) v = l + e(p 7 ) . 
Since cycle point (I, u) lies in A 8 , applying Lemma [3741 

(170) u = I + e(p 8 ) . 

From (11631 ) and (11651) eq. 1 17 H is obtained and substituting the other equations into 
this the desired result is found. 

(171) x + e{p l ) = t + e(p 3 ). 

□ 

It has hence been proved that the conditions listed on the statement considers 
all the possible 8-cycles that can exist on the studied quasi-cyclic matrices. 

□ 
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